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1  INTRODUCTION 


The  target  resolution,  parameter  estimation  and  clutter  rejection  capabilities  of 
a  radar  are  determined  by  the  ambiguity  fimction  associated  with  the  transmitted 
waveform  [1-3].  The  ambiguity  function  and  its  plot  provide  a  qualitative  assessment 
of  the  suitability  of  different  waveforms  in  meeting  system  requirements.  Theory 
and  application  of  the  ambiguity  function  to  monostatic  radars  in  a  Gaussian  noise 
environment  is  fairly  well  understood  [4].  A  monostatic  radar  refers  to  a  radar  system 
where  the  transmitter  and  the  receiver  are  located  at  the  same  site.  Over  the  past 
several  years  there  has  been  quite  a  bit  of  interest  in  the  area  of  bistatic  and 
multistatic  radars.  This  is  primarily  due  to  survivability  considerations.  In  these 
systems,  the  transmitters  and  the  receivers  are  located  at  different  sites  in  a  variety  of 
configurations.  The  concept  and  application  of  ambiguity  function  for  bistatic  radars 
has  not  been  examined  in  detail.  One  of  the  main  objectives  of  this  study  was  to 
investigate  the  concept  of  ambiguity  function  in  the  context  of  bistatic  radars.  Before 
discussing  the  bistatic  radar,  it  is  useful  to  review  the  manner  in  which  the  ambiguity 
fimction  arises  in  a  monostatic  radar.  In  Chapter  2,  we  present  the  derivation  of 
ambiguity  fimction  from  the  basic  principles  of  detection  and  estimation  theory.  In 
Chapters  3  and  4,  we  present  the  formulation  of  ambiguity  function  for  bistatic  radars. 
In  Chapter  3,  we  describe  the  North-referenced  coordinate  system  for  the 
representation  of  bistatic  geometries.  Procedures  for  transformation  from  the  Global 
Earth-centered  coordinate  system  and  the  local  north-east-down  coordinate  system  to 
the  North-referenced  coordinate  system  are  also  given.  This  will  enable  the  use  of  our 
methodology  even  when  the  measurements  are  obtained  in  these  other  coordinate 
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systems.  In  Chapter  4,  we  discuss  in  detail  the  effect  of  bistatic  geometry  on 
ambiguity  function  analysis.  Two  examples  are  provided  to  illustrate  the  nature  of 
ambiguity  fimction  for  a  variety  of  system  configurations.  These  examples  consider 
a  Gaussian  pulse  and  a  waveform  consisting  of  three  rectangular  pulses. 

Most  of  the  development  of  ambiguity  function  is  based  on  the  processing  of 
continuous  waveforms  at  the  receiver.  However,  modern  radars  employ  digital  signal 
processing  techniques  in  that  the  received  signals  are  sampled  prior  to  signal 
processing.  It  is,  therefore,  imperative  that  the  ambiguity  function  in  terms  of 
discrete-time  signals  be  analyzed  in  detail.  In  Chapters  5-8,  we  examine  a  number  of 
issues  related  to  the  ambiguity  function  based  on  sampled  signals.  In  Chapters  5  and 
6  we  derive  expressions  for  the  continuous-time  continuous-frequency  ambiguity 
fimction  based  on  sampled  signals.  In  Chapter  7,  the  relationship  between  the 
discrete-time  discrete-frequency  and  continuous-time  continuous-frequency  ambiguity 
fimctions  is  established.  Sampling  rates  required  to  avoid  aliasing  are  also 
determined.  In  Chapter  8,  the  reconstruction  of  the  conventional  ambiguity  function 
from  the  discrete-time  discrete-frequency  ambiguity  fimction  is  discussed. 

One  of  the  important  problems  in  the  operation  of  a  bistatic  radar  is  that  of 
synchronization,  i.e.,  the  availability  of  the  transmitted  waveform  replica  at  the 
receiver.  In  the  cooperative  mode  of  operation,  this  availability  is  assumed.  In 
Chapter  4,  ambiguity  function  analysis  for  a  cooperative  bistatic  radar  was  carried  out. 
When  the  bistatic  radar  operates  in  a  noncooperative  mode,  uncertainties  exist  in  the 
reference  signal  available  at  the  receiver.  In  Chapter  9,  this  problem  is  briefly 
considered  and  a  signal  model  that  includes  these  uncertainties  is  suggested.  In 
Chapter  10,  a  summary  of  the  major  contributions  of  this  investigation  is  provided. 
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2  AMBIGUITY  FUNCTION  FOR  MONOSTATIC  RADARS 


This  chapter  is  devoted  to  the  study  and  derivation  of  ambiguity  functions  for 
monostatic  radars.  It  is  shown  that  the  notion  of  ambiguity  function  comes  up 
naturally  when  dealing  with  the  detection  and  parameter  estimation  problems.  It  is 
typically  derived  assuming  a  slowly  fluctuating  point  target.  The  model  associated 
with  such  a  target  is  developed  next 

2.1  Slowly  fluctuating  point  target  model 

In  the  propagation  direction  of  an  electromagnetic  pulse,  assume  a  target  has 
scatterers  distributed  over  the  range  interval  (Rj,  R2),  as  shown  in  Fig.  2-1.  Let  tj 
denote  the  time  instant  at  which  the  leading  edge  of  the  pulse  reaches  range  Rj  where 
i  =  1,2.  We  assume  the  time  origin  to  correspond  to  the  instant  at  which  the  leading 
edge  of  the  pulse  appears  at  the  transmitting  antenna.  If  T  denotes  the  time  duration 
of  the  transmitted  pidse,  we  say  that  the  target  can  be  modeled  as  a  point  target 
provided 

tj  -  tj  <  <  T. 

Assiune  the  transmitted  pulse  contains  n  complete  cycles  of  carrier.  (Here, 
"transmitted  pulse"  refers  to  a  single  finite-duration  radar  waveform  which  may 
actually  be  composed  of  a  number  of  simple  subpulses.)  It  follows  that  the  carrier 
frequency  is  related  to  the  pulse  width  by 


Since 
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=  c, 

where  X  and  c  are  the  wavelength  and  the  velocity  of  propagation,  respectively,  of  the 
transmitted  pulse,  it  follows  that 


X  =  —  =  — 


At  time  t,  assume  that  the  leading  edge  of  the  transmitted  pulse  arrives  at  range  R. 


Since 


R  =  ct, 


the  time  difference  (t^-ti)  can  be  expressed  as 


h 


R,-R, 


Thus,  the  inequality  used  to  define  a  point  target  becomes 


We  conclude  that  a  target  can  be  modeled  as  a  point  target,  that  is,  the  return  pulse 
shape  closely  approximates  the  transmitted  pulse,  provided 

R^-Rj  <  <  nX. 

It  should  be  noted  that,  if  the  transmitted  pulse  consists  of  a  sequence  of  several 
subpulses,  n  still  refers  to  the  number  of  carrier  cycles  in  one  pulse.  For  convenience, 
some  people  use  the  more  restrictive  condition 

Rj  -  Rj  <  <  ^ 


in  order  to  define  a  point  target. 
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To  examine  the  significance  of  the  tenn  "slowly  fluctuating", 

transmitted  pulse  is  given  by 

assume  that  the 

8^.(0  =  sjl  Aft)  COS[0)^t  +  <|)(t)] 

=  v/2  Re{^  A(t)e^‘)e^“‘=‘} 

=  y/l  Re{y^  f(t)eJ““‘) 

where  Re  {•}  denotes  the  real  part  operation. 

(2.10) 

o>  =  2jif , 

C  c* 

(2.11) 

and 

f(t)  =  A(t)e^‘> 

(2.12) 

is  the  complex  envelope  of  the  transmitted  pulse.  Let  the  complex  envelope  be 

normalized  such  that 

/“  |f(t)|2  dt  =  1. 

(2.13) 

Then  the  energy  of  the  transmitted  signal  is  equal  to  E,. 

For  a  monostatic  radar,  assume  that  a  point  target  is  located  at  range  with 

radial  velocity  v,.  The  geometry  is  indicated  in  Fig.  2-2.  Let  the  point  target  consist 

of  a  rough  surface  with  K  scatterers.  The  received  waveform  due  to  the  reflection 

from  the  i“'  scatterer  is 

Sj  (t)  =  V2  Re  {  ^  } 

(2.14) 

where  gj  and  0,  denote  the  attenuation  and  phase  shift,  respectively,  incurred  in  the 

reflection  process,  x,  =  2Ra/c  denotes  the  round-trip  time  delay  and 

6 

“d.  =  —  “c 


(2.15) 


is  the  shift  in  the  carrier  frequency  due  to  the  radial  velocity  of  the  target.  cOoa  is 
called  the  Doppler  shift.  Summing  the  returns  from  the  K  scatterers,  the  total 
received  waveform  is 


SR(t)  =  Y,  s,(t) 

1=1 

=  ^2  Re  I  yi;  £  g,  e«  f  (t- ^ 

=  ^  He  { 6  f 

where  the  operations  of  summation  and  real  part  have  been  interchanged  and 


(2.16) 


K  k 

^  =  E  8ie^‘=E«i' 

i=l  i=l 


(2.17) 


Assiuning  that  the  complex  random  numbers  g.  are  statistically  independent  and  that 
K  is  sufficiently  large,  we  can  use  a  central  limit  theorem  argument  to  conclude  that 


6  =  bi  +  jbg  = 

is  a  zero-mean  complex  Gaussian  random  variable.  It  follows  that 


(2.18) 


|6|  =  +  bg 


(2.19) 


is  a  Rayleigh  distributed  random  variable  and 


05  =  tan 


-1 


(2.20) 


is  a  imiformly  distributed  random  variable.  Absorbing  the  constant  factors 
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into  the  phase  component  of  b,  the  total  received  signal  becomes 


s,(t)  .  J2  Ref/i;  6  f(t-Ty"”-‘ 


(2.22) 


The  point  target  is  defined  to  be  slowly  fluctuating  if  b  can  be  asstuned  to  be 

constant  while  the  target  is  illuminated  by  the  transmitted  pulse  of  duration  T. 

Having  defined  what  is  meant  by  a  slowly  fluctuating  point  target,  we  now 
show  how  the  ambiguity  function  arises  in  the  radar  detection  problem. 

2.2  The  Radar  Detection  Problem 

Let  f(t)  denote  the  complex  envelope  of  the  transmitted  radar  pulse.  Assume 

that  f(t)  is  zero  outside  the  interval  0  <  t  <  T.  Given  the  complex  envelope  of  the 

received  waveform  over  the  interval  Xh  <  t  <  %  +  T,  we  are  interested  in  determining 
whether  a  slowly  fluctuating  point  target,  at  range  R  with  radial  velocity  v,  is  located 
at  the  point 


,  2R  2v  . 

(t„  =  —  , 


(2.23) 


in  the  delay  -  Doppler  plane.  In  the  absence  of  a  target,  the  complex  envelope  of  the 
received  waveform  is  assumed  to  be 


f(t)  =  n(t),  Tjj  t  +  T  (2-24) 

where  n(t)  denotes  the  complex  envelope  of  a  Gaussian  white  noise  process.  On  the 
other  hand,  if  a  slowly  fluctuating  point  target  is  located  at  the  point  (iH,cOnJ  in  the 


delay  -  Doppler  plane,  the  complex  envelope  of  the  received  waveform  is 

f(t)  =  ^  6  f(t-Tjj)e^“”«‘  +  n(t),  Tjji  t  i  Tjj  +  T  .  (2.25) 

Hence,  given  the  two  hypotheses 


Hq  :  f(t)  =  n(t) 


Hj  :  f(t)  =  ^  6  f(t-tj,)e^“”«‘+n(t), 


T„itiT„+T 


(2.26) 


we  wish  to  determine  the  Neyman-Pearson  receiver,  i.e.,  the  receiver  that  maximizes 
the  detection  probability  for  a  specified  false  alarm  probability,  for  detecting  a  slowly 
fluctuating  point  target  at  the  hypothesized  point  (Xh,  cOd)  in  the  delay-Doppler  plane. 

For  this  purpose,  it  is  convenient  to  expand  f(t)  in  a  Karhunen-Loeve  expansion 
as  given  by 


r(0  =  X)  ^  t  +  T.  (2.27) 

i=l 

Because  b  is  a  complex  Gaussian  random  variable  and  n(t)  is  a  white  complex 
Gaussian  random  process,  the  coefficients 

fj  =  m  4>;(t)  dt,  (2.28) 

are  statistically  independent  complex  Gaussian  random  variables.  In  addition,  since 
the  noise  is  white,  the  basis  functions  0.(t)  can  be  chosen  from  any  set  which  is 

complete  and  orthonormal  over  the  interval  Xh  ^  t  <  Xh  +  T.  For  convenience,  let  the 
first  basis  function  be 
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2.3  Derivation  Based  on  Sufficient  Statistics 

The  remaining  orthonormal  functions  in  the  basis  set  can  be  selected  in  any 

convenient  manner.  For  i  >  1,  the  ith  coefficient  in  the  Karhunen-Loeve  expansion  of  f(t) 
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is 


(2.35) 


where  fi;  is  the  ith  coefficient  in  the  Karhunen-Loeve  expansion  of  n(t).  Since  n(t)  is 
a  complex  Gaussian  white  noise  process,  fl;  for  i  >  1  is  statistically  independent  of  nj. 


Also,  by  assumption,  hj  is  statistically  independent  of  b.  Thus,  f.  for  i  >  1  contains 


no  information  either  about  the  target  return  or  the  noise  coefficient  fii.  It  follows 
that  fj  is  a  sufficient  statistic.  The  log  likelihood  ratio  is  given  by 


h  .  -3^|R.P  .  (2.36) 

P,,|H.(«imj)J  [2ojE,tN„j  No(2o^,+N^ 

Because  the  constants  can  be  combined  with  the  threshold,  the  Neyman-Pearson 
receiver  performs  the  likelihood  ratio  test 


H, 

|R.P  :  y  (2.37) 

Ho 

where  the  threshold  y  is  selected  so  as  to  achieve  a  specified  false  alarm  probability. 

The  above  test  maximizes  the  detection  probability  under  a  false  alarm 
probability  constraint  assuming  that  a  slowly  fluctuating  point  target,  if  present,  is 
located  at  the  point  (Xh/  in  the  delay  -  Doppler  plane.  However,  the  target  may 

actually  be  located  at  the  point  (x^  ,  cou,).  The  complex  envelope  of  the  received 
waveform  is  then 
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f(t)  =  ^  6  n(t). 


(2.38) 


Hence, 


Jtn 


(2.39) 


It  follows  that 


|fip  =  Ejb|2|/^"^''f(t-T,)f‘(t-XH)e^““‘‘““«^dtp 
+  tmns  involving  n(t). 


(2.40) 


Ignoring  the  noise  terms  and  the  multiplicative  factor  EJb|^,  the  ambiguity  fimction 


is  defined  to  be 


(2.41) 


Because  f(t-x„)  is  zero  outside  the  interval  %  <  t  <  Xh  +  T,  the  range  of  integration 

H 

may  be  extended  from  -  oo  to  +  «».  Hence,  the  ambiguity  function  is  commonly 
expressed  as 


(2.42) 


It  is  convenient  to  express  the  ambiguity  function  in  terms  of  the  differences 
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x'  =  and  o>4  = 


This  results  in 


Moving  the  factor 


(2.43) 


=  l/’ f(t-Tj,-T')r(t-Tjj)e^““^‘‘"«*"“^dtp.  (2.44) 


J(«'d)Th 


(2.45) 


outside  of  the  integral  and  noting  that 


we  have 


(2.46) 


0(Th,  t,,  QdP  =  ||_*  f(t-Tg-T')f’(t-TH)e^““^‘’""^dtp 


(2.47) 


A  shift  of  the  time  origin  results  if  t-Xn  is  replaced  by  L  This  gives 


0(Th,x,,(Od^,Wd,)  =  “d) 

=  |/“  fft-xOfW'^dtp. 


(2.48) 


A  symmetrical  form  of  the  ambiguity  fxmction  is  obtained  by  a  further  shift  of  the 
time  origin.  Replacing  t  by  t  +  x72  gives 


0(x'.<J^,)  =  l/j(t-y)fVy)  e  2  Jit|2 


(2.49) 


Moving  the  factor  e^“  outside  of  the  integral  and  noting  that 


14 


(2.50) 


the  symmetrical  form  of  the  ambiguity  function  reduces  to 


(2.51) 


The  effect  of  changing  over  to  variables  x'  and  (%'  is  to  shift  the  target  to  the  origin 


of  the  x',cDd'  plane.  In  other  words,  when  x'  =  (Op'  =  0,  then 


Th  =  r.  and 


(2.52) 


Clearly,  to  minimize  the  likelihood  that  a  target  will  be  detected  at  the  point  Xjj,C0p^) 

when  it  is  actually  located  at  the  point  (x^,  cou,),  it  is  desirable  that  0(x',  co  p)  be  an 
impulse  located  at  the  origin.  This  ideal  ambiguity  function  cannot  be  achieved  in 
practice. 

2.4  Alternative  Derivation 

It  is  instructive  to  rederive  the  ambiguity  function  by  an  alternative  approach 
to  the  detection  problem.  In  the  previous  development,  use  was  made  of  the  fact  that 
the  first  coefficient  in  the  Karhunen-Loeve  expansion  of  the  complex  envelope  of  the 
received  signal  is  a  sufficient  statistic.  Even  so,  the  solution  could  have  proceeded  by 
considering  all  of  the  coefficients  in  the  expansion 


f(t)  =  E  ■'h  ^  t  ^  +T. 

i=l 


(2.53) 


This  latter  approach  is  discussed  next 

As  before,  the  two  hypotheses  are  assumed  to  be 
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(2.54) 


Hq  :  f(t)  =n(t) 

Hj  ;  f(t)  =  s(t)  +  n(t)  ^  t  ^  ^  T 


where 


s(t)  =  ^  6  f(t-TH)e^““«‘ 


(2.55) 


Recalling  that  E[  I  b  h  ]  =  2a^b  /  the  covariance  function  of  s(t)  is  given  by 


K,(t,u)  =  E[s(t)s*(u)I  =  f(t-TH)r(u-tH)e 


(2.56) 


where 


i,  =  2ajE.  . 


(2.57) 


The  complex  envelope  of  the  Gaussian  white  noise  process  is  assumed  to  be 
statistically  independent  of  s(t)  with  covariance  function 


K,(t,u)  =  E[n(t)n*(u)]  =  N„6(t-u). 


(2.58) 


Keeping  all  of  the  terms  in  the  Karhunen-Loeve  expansion  of  f(t),  the  sufficient 
statistic  for  the  Neyman-Pearson  receiver  is  given  by  [4] 


«  =  f*(t)h(t,u)f(u)dtdu 


(2.59) 


where  h(t,u)  is  the  solution  to  the  integral  equation 


Noh(t,u)  +  j  h(t;z)Kg(z,u)dz  =  Kj(t,u),  Tjj^t,uiTjj+T. 


(2.60) 


Substituting  for  the  signal  covariance  fimction,  the  integral  equation  becomes 
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(2.61) 


Noh(t,u)  +  f  h(t^)E^f(z-Tjj)f*(u-TH)e^“°®^^  ”^dz 


=  E^f(t-XH)f’(u-TH)e 


By  inspection,  the  solution  is  seen  to  be 


h(t,u) 


;,u)  =  ^3—  f(t-Tjj)r(u-TH)e^“®"‘‘ 


(2.62) 


where  use  is  made  of  the  fact  that  f(t)  is  normalized  to  have  unit  energy  such  that 


/  =  E 


(2.63) 


Substitution  of  the  expression  for  h(t,u)  into  the  sufficient  statistic  results  in 


{  =  J_  ff  f*(t)  fCt-Tn)  f ’(u-TH)e^'^^‘ “V(u)dtdu 

N„ 


(2.64) 


- 3 -  /  f(t)r(t-Xjj)e 

No(E,-No)  i 


No(E,+No) 

where  is  the  first  coefficient  in  the  Karhunen-Loeve  expansion  of  f(t)  as  defined 

previously.  Because  the  multiplicative  constant  can  be  combined  with  the  threshold, 
the  above  result  is  consistent  with  the  likelihood  ratio  test  derived  in  the  previous 
approach.  The  second  approach  was  discussed  because  it  proves  to  be  more  useful 
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when  considering  sampled  data. 

2.5  The  Parameter  Estimation  Problem 


The  ambiguity  function  is  also  related  to  the  parameter  estimation  problem  for 
a  slowly  fluctuating  point  target  where  the  complex  envelope  of  the  received  signal 
is  assumed  to  be  of  the  form 


f(t)  =  6  Tjj  i  t  i  Tjj  +  T 


(2.65) 


and  Xh  and  are  unknown  nonrandom  parameters  that  we  wish  to  estimate. 

Because  the  unknown  parameters  are  nonrandom,  we  use  a  maximum 
likelihood  estimation  procedure.  The  maximum  likelihood  estimates  are  those  values 

of  Xh  and  for  which  the  likelihood  function  is  a  maximiun.  The  derivation  of  the 

likelihood  function  is  similar  to  that  of  the  likelihood  ratio.  The  first  step  is  to  make 
a  Karhunen-Loeve  expansion  for  f(t)  assuming  Xjj  and  are  known.  This  means 

that  we  must  choose  basis  fimctions  that  are  conditionally  dependent  on  Xh  and  . 


We  write 


m  =  E  ri(TH.WDH)  ^  ^ 


(2.66) 


where 


f(t)^;(t:xj,,a)p^  dt 


(2.67) 


For  convenience,  let  the  first  basis  fimction  be 
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(2.68) 


The  first  coefficient  in  the  expansion  of  f(t)  is  then 


(2.69) 


where 


(2.70) 


is  the  first  coefficient  in  the  Karhunen-Loeve  expansion  of  the  complex  white  noise 
process  assuming  Th  and  are  given.  Because  b  is  a  complex  Gaussian  random 

variable  and  n(t)  is  a  complex  Gaussian  random  process,  the  coefficients  f.  are 
statistically  independent  complex  Gaussian  random  variables.  As  before,  let 


E[|nip]  =  and  E[|6|2]  =  20^. 
The  conditional  probability  density  function  of 


(2.71) 


is  then  given  by 


(2.72) 


J _  g  2oJe,*N„ 


7i(2oJe,+No) 


(2.73) 


The  remaining  orthonormal  fimctions  in  the  basis  set  can  be  selected  in  any 
convenient  manner.  For  i  >  1,  the  ith  coefficient  in  the  Karhunen-Loeve  expansion  of  f  (t) 
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IS 


where  ni(xH,Ci>jj^)  is  the  ith  coefficient  in  the  Karhunen-Loeve  expansion  of  n(t) 
asstuning  Xh  and  Wp  are  known.  Because  the  coefficients  fiiixn,  Wj,  ),  for  i  >  1,  depend 

H  H 

only  on  the  noise  and  are  statistically  independent  of  njlxH/Wjj^),  they  contain  no 
useful  information  for  estimating  Xh  and  We  conclude  that  fj(Xj^,C0p  )  is  a 

sufficient  statistic.  Since  the  conditional  probability  density  function 

I  <2-75) 

is  maximized  by  maximimizing  the  exponent,  it  follows  that  the  maximum  likelihood 
estimates  of  Xh  and  Uj,  are  those  values  of  Xh  and  Wp  for  which 

H  H 

m  (2.76) 

is  a  maximum. 

As  in  the  detection  problem,  the  ambiguity  fimction  is  defined  by  ignoring 
those  terms  in  |Rj(Xj^,C0p  )|^  that  involve  n(t).  If  the  actual  values  of  delay  and 

Doppler  frequency  are  x^  and  Up  ,  the  ambiguity  function  is  found  to  be 


=  l/^  f(t-x,)r(t-Tjj)e^^“°‘ (2.77) 


which  is  identical  to  that  obtained  in  the  detection  problem.  Thus,  introducing  the 


20 


change  of  variables 

z'  =  T,-Tjj  and  co'p  = 

the  ambiguity  function  can  again  be  written  in  the  symmetrical  form 

=  |/“  f{t-— S  f*(t+— )  e^“®‘dtp.  (2-79) 

Clearly,  to  minimize  the  maximum  likelihood  estimation  errors,  it  is  desirable  that 
0(x',(o'd)  be  an  impulse  located  at  the  origin.  As  pointed  out  previously,  this  ideal 
ambiguity  function  cannot  be  achieved  in  practice. 

As  opposed  to  simply  defining  the  ambiguity  function,  as  is  done  in  most 
references,  we  have  shown  that  the  ambiguity  function  arises  naturally  in  both  the 
detection  and  parameter  estimation  problems.  In  both  applications,  it  is  desirable  to 

choose  a  waveform  for  f(t)  such  that  the  shape  of  the  ambiguity  fimction  in  its 

symmetrical  form  approximates  an  impulse  located  at  the  origin.  Examples  and 
properties  of  the  ambiguity  fimction  are  well  documented  and  are  not  repeated  here. 
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3  COORDINATE  SYSTEMS 


Much  of  the  radar  systems  analysis  employs  the  Global,  Earth-centered 
coordinate  system  using  the  earth  as  the  reference.  This  three-dimensional  coordinate 
system  was  also  employed  in  the  previous  Rome  Laboratory  efforts  on  bistatic  radar. 
However,  this  coordinate  system  is  quite  cumbersome  and  unnecessarily  complex  for 
expressing  the  bistatic  configxiration.  Jackson  [5]  presented  the  North-referenced 
coordinate  system  to  describe  the  bistadc  radar  geometry.  It  is  a  two-dimensional 
system  and  is  fairly  convenient  to  use.  Willis  [6]  further  elaborated  on  this  system 
and  has  employed  it  extensively  to  analyze  the  bistatic  radar  system  and  to  obtain 
many  useful  results.  For  this  reason,  we  also  employ  the  North-referenced  coordinate 
system  for  the  ambiguity  function  analysis  for  the  bistatic  radar.  In  this  section,  both 
of  these  coordinate  systems  are  introduced.  The  procedure  for  transformation  from 
the  Global,  Earth-centered  coordinate  system  to  the  North-referenced  coordinate 
system  is  derived.  This  will  enable  the  use  of  our  results  if  the  measurements  are 
provided  in  the  Global,  Earth-centered  coordinate  system. 

Another  common  measurement  practice  is  to  use  a  local  coordinate  system.  An 
example  of  such  a  coordinate  system  is  the  north-east-down  coordinate  system.  We 
present  this  coordinate  system  and  the  procediue  for  converting  the  measurements 
provided  in  this  system  to  the  North-referenced  system.  If  any  other  coordinate 
systems  are  employed  for  measurements,  appropriate  coordinate  conversion 
procedures  can  be  developed  so  that  our  results  can  be  used. 

3.1  North-referenced  Coordinate  System 

This  is  a  two-dimensional  coordinate  system  confined  to  the  bistatic  plane 
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Tx  -  Transmitter 
Rx  -  Receiver 
Tgt  -  Target 
N  -  North  of  North- 

referenced  Coordinate 

System 


9j  =  transmitter’s  look  angle 

6r  =  receiver’s  look  angle 

(3  =  8x-0r 

L  =  baseline  range 

Rx  =  transmitter  to  target  range 

Rr  =  receiver  to  target  range 

R  =  Rx+Rr 


Fig.  3-1:  North-referenced  Coordinate  System 
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formed  by  the  transmitter  the  receiver  R,  and  the  target  Tgt.  Figure  3-1  shows  the 
coordinate  system  and  its  parameters.  This  coordinate  system  is  represented  as  in 
maps  with  a  north-up  orientation.  The  distance  between  the  transmitter  and  the 
receiver  is  the  baseline  range,  or  simply  baseline  and  is  denoted  by  L.  The  receiver 
is  at  the  eastern  end  and  the  transmitter  is  at  the  western  end  of  the  baseline.  The 
triangle  formed  by  the  transmitter,  the  receiver  and  the  target  is  called  the  bistatic 
triangle.  Transmitter's  look  angle  0t  and  receiver's  look  angle  6*  are  measured  from 
the  north  of  the  coordinate  system,  positive  clockwise.  Also,  the  angles  Sj  and  6*  are 

restricted  to  the  interval  —1.  The  bistadc  angle  p  is  the  angle  between  the 

2  2 

transmitter  and  the  receiver  in  the  bistatic  triangle  with  the  target  as  the  vertex.  It  is 
given  by  p  =  0x  •  0R.  Using  the  definition  and  the  restriction  on  the  values  Sj  and  6r, 
the  bistatic  angle  is  positive  and  in  the  interval  [0,  k]  when  the  target  is  to  the  north 
of  the  baseline,  and  is  negative  and  in  the  interval  [-n,  0]  when  the  target  is  to  the 
south  of  the  baseline.  The  line  that  bisects  the  bistatic  angle  is  called  the  bistatic 
bisector.  The  region  to  the  north  of  the  baseline  is  called  the  northern  hemisphere, 
and  likewise,  the  southern  hemisphere  is  the  region  to  the  south  of  the  baseline.  In 
general,  bistatic  radar  operation  and  performance  in  the  northern  and  southern 
hemispheres  are  equivalent  for  symmetric  geometries.  The  range  from  the  transmitter 
to  the  target  is  denoted  by  Rj.  The  range  from  the  receiver  to  the  target  is  denoted 
by  Rr.  For  a  bistatic  radar,  the  target  range  is  measured  in  terms  of  the  total  range, 
or  range  sum,  given  by  R  =  Rx  +  Rr.  Note  that  in  bistatic  radar  a  constant  range 
contour  can  be  described  by  an  ellipse.  This  is  defined  as  the  locus  of  points  for 
which  Rx  +  Rr  is  constant  For  a  monostatic  radar,  the  constant  range  contour  reduces 
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to  a  circle  because  in  this  case,  Rt  =  Rr.  An  example  of  constant  range  or  constant 
delay  contours  is  shown  in  Fig.  3-2.  Also,  constant  signal-to-noise  ratio  contours  can 
be  described  in  terms  of  ovals  of  Cassini  These  ovals  are  the  loci  of  points  for  which 
the  product  Rt  Rr  is  constant  Typical  ovals  of  Cassini  are  shown  in  Fig.  3-3.  Again, 
for  a  monostatic  radar,  the  constant  SNR  contour  reduces  to  a  circle.  Note  that  in 
bistatic  radars,  these  ovals  of  Cassini  define  three  distinct  operating  regions:  (I) 
receiver-centered  region,  the  small  oval  around  the  receiver  as  given  in  Fig.  3-3;  (II) 
transmitter-centered  region,  the  small  oval  arround  the  transmitter;  and  (III) 
receiver/transmitter-centered  region,  called  the  cosite  region,  any  of  the  regions 
surroimding  both  the  transmitter  and  the  receiver.  As  given  in  [6],  the  receiver 
centered  region  can  be  used  for  applications  such  as  passive  situation  awareness 
where  one  is  interested  in  the  detection  of  a  target  around  a  passive  receiver.  The 
transmitter-centered  region  can  be  used  with  a  non-cooperative  transmitter.  The  cosite 
region  can  be  used  for  medium  range  air  defense  and  satellite  or  missile  tracking. 

The  bistatic  radar  geometry  can  be  defined  in  terms  of  foxir  parameters,  namely, 
0T,  0R,  L  and  R.  However,  it  is  sufficient  to  know  any  three  of  these,  and  the  fourth 
one  can  be  obtained  from  the  knowledge  about  the  other  three.  Also,  Rt  and  Rr  can 
be  determined.  In  the  following,  relations  among  the  parameters  are  developed. 
Three  cases  are  considered. 

Case  Ii  0n.  R  and  L  are  specified. 

In  this  case,  we  assume  that  the  total  range  R,  baseline  L  and  the  receiver's  look 
angle,  or  the  angle  of  arrival  of  the  target  return,  0r,  are  available.  Applying  the  law 
of  cosines  to  the  bistatic  triangle,  we  have 
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Fig.  3-2:  Concentric  Ellipses  (Constant  Range  Contours) 


26 


operating  regions:  I  transmitter-centered  region 
n  receiver-centered  region 
111  cosite  region 


Fig.  3-3:  Contours  of  Constant  Signal-to-Noise  Ratios  (Ovals  of  Cassini) 


rJ  =  Ri  +  -  2R^L  cos  (90°  +  6^) 

By  definition,  Rj  =  R  -  Rr.  Then, 

(3.1) 

Rt'  =  (R  -  Rr)*  =  R"  -  2RRr  +  Rr* 

Combining  (3.1)  and  (3.2),  we  have 

(3.2) 

R"  -  2RRr  +  Rr"  =  Rr'  +  L'  -  2Rr  L  cos  (90®+eR) 

Reananging  terms  of  the  two  sides  of  (3.3),  we  obtain 

(3.3) 

p  - 

^  2(R+Lsin6R) 

From  (3.1), 

(3.4) 

Rt  =  [Rr^  +  L*  +  2RRLsineR]'' 

Using  (3.4),  we  may  obtain  another  expression  for  Rj 

(3.5) 

R2+L2+2RLsineR 

K.  =  R-R,  = - — 

^  ^  2(R+Lsm0R) 

(3.6) 

Now,  the  distance  between  the  target  and  both  the  transmitter  and  the  receiver  have 

been  obtained. 

Again,  using  the  law  of  cosines,  it  can  be  written  that 

Rg^  =  R/  +  L*  -  2RxL  cos  (90®  -  Gx) 

=  (R  -  Rr)"  +  L*  -  2(R  -  RR)Lsinex 

=  R"  -  2RRr  +  Rr*  +  L*  -  2(R  -  RR)Lsinex 

Rearranging  the  two  sides  of  (3.7),  we  have 
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(3.7) 

(3.8) 


sin6^ 


r2+l2-2RR^ 

2(R-R^)L 


Substituting  (3.4)  into  (3.8),  we  have 


.  ,  (R*»LVn0.*2aL 

— — — - , 


7c-sin 


(Jl^+L2)8in0j+2IIL 

(R*+L^)+2RL8iii0a’ 


(3.9a) 


(3.9b) 


where  --  i  sin'*(0  <  -•  Note  that  (3.9a)  represents  the  case  when  the  target  is  in  the 
2  ''^2 

northern  hemisphere,  i.e.,  -  |  ^  6*  <  The  operation  sinM •)  yields  0t  which  is  also 

in  the  northern  hemisphere.  However,  when  the  target  is  in  the  southern  hemisphere, 
Le.,  -  i  0R  <  y.  0T  “  defined  as  in  (3.9b)  so  that  the  resulting  0t  also  corresponds 

to  the  southern  hemisphere. 

Thus,  using  (3.4),  (3.6)  and  (3.9),  the  bistatic  triangle  is  specified  and  the  target's 
position  is  determined  with  the  measurements  of  0ii,  R  and  L. 

Case  IT:  0t.  R  and  L  are  specified 

In  this  case,  the  total  range  R,  baseline  L  and  the  transmitter's  look  angle  0t  are 
assumed  to  be  given.  Applying  the  law  of  cosines  to  the  bistatic  triangle,  we  have 

Rr"  =  Rt*  +  L"  -  2RtL  cos  (90“  -  0t)  (3.10) 

By  definition  Rr  =  R  -  Rj.  Then 

Rr*  =  (R  -  Rt)*  =  R*  -  2RRt  +  Rt* 
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Combining  (3.10)  and  (3.11), 

R*  -  2RRt  +  Rt*  =  Rj  +  V  -  2Rx  L  cos  (90"  -  Qj) 
Rearranging  terms  of  the  two  sides  of  (3.12),  we  obtain 


(3.12) 


R"-L" 

2(R-Lan0-) 


(3.13) 


From  (3.10), 


R.  =  [R  *  +  L*  .  2RxLsineT]'' 


(3.14) 


Alternatively,  using  (3.13),  we  may  obtain  another  expression  for  Rg, 


Rj  =  R  -  Rj. 


R2+L2-2RLsme, 

2(R-Lsm6T) 


(3.15) 


Thus,  the  distances  between  the  target  and  both  the  transmitter  and  the  receiver  have 
been  obtained. 

Next,  using  the  law  of  cosines,  it  is  also  found  that 
R/  =  Rg2  +  L*  -  2RrL  cos  (90"  +  Or) 


=  (R  -  Rt)"  +  L'  -  2(R  -  Rt)L  cos  (90"  +  0r) 
=  R*  -  2RRt  +  Rt*  +  L*  +  2(R  -  Rt)  LsinOR 
Rearranging  the  two  sides  of  (3.16),  we  have 


(3.16) 


sme^  = 


2RRr-R2-L= 

2(R-Rr)L 


(3.17) 


Substituting  (3.13)  into  (3.17),  we  have 
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(3.18a) 


,  (R2+L2)sm0T.-2RL 

sin  — ; — ; - . 

R^+L^-2RLsine~ 


--  i  0T  <  - 
2^2 


,  (R2+L2)sm0T.-2RL 

Tc-sin"*  — - — - , 

R2+L2-2RLain0^ 


-  s  0T  <  ^ 
2^2 


(3.18b) 


As  before,  the  cases  of  northern  and  southern  hemispheres  are  treated  separately  in 
(3.18a)  and  (3.18b)  respectively. 

Hence,  using  (3.13),  (3.15)  and  (3.18),  the  bistatic  triangle  is  specified  and 
target's  position  is  determined  with  the  measurements  6x,  R  and  L. 


In  this  case,  the  angle  measurements  Sr  and  &r,  and  the  range  measurement  R 
are  assumed  to  be  available.  Consider  Fig.  3*4.  It  is  found  that,  in  the  bistatic 


triangle. 


Rj  cos  6t  =  Rr  cos  6r 


Substituting  Rx  =  R  -  Rr  into  (3.19), 


(R  -  Rr)  cos  0t  =  Rr  cos  0r 


Rearranging  the  two  sides  of  (3.20),  we  have 


R  cos  0x 
cos  0K^<^ 


(3.19) 


(3.20) 


(3.21) 


Rp  =  R  Rn  ~ 


R  cos  0R 
cos  0|^-t-cos  Op 


(3.22) 
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Fig.  3-4:  Bistatic  triangle  and  associated  parameters 


(3.23) 


In  the  bistatic  triangle,  it  is  also  true  that 

L  =  Rt  sin  St  -  Rr  sin  6* 

Substituting  (3.21)  and  (3.22)  for  Rr  and  Rt,  respectively,  into  (3.23),  we  have 

^  ^  RsinCe-r-e^)  0.24) 

cos  6^-t-cos8|i 


Now,  with  (3.21),  (3.22)  and  (3.24),  and  the  given  angle  measurements,  the  bistatic 
triangle  is  totally  specified. 

3  2  Farfli-ggntered  Coord»"afa»  System  (POD  Worid  Geodetic  Systeml 

This  is  a  three-dimensional,  right-handed  rectangular  coordinate  system,  with 
its  origin  set  at  the  center  of  the  earth,  as  shown  in  Figure  3-5.  The  z-axis  pomts  to 
the  true  north  of  the  earth,  and  the  x-y  plane  contains  the  equator.  The  x-axis  passes 
through  the  zero  (Greenwich)  meridian.  Position  of  an  object  is  often  provided  in  this 
coordinate  system.  Sometimes  it  is  more  convenient  to  specify  the  position  of  an 
object  in  terms  of  its  geographic  coordinates,  Le.,  the  position  is  given  in  terms  of  its 
latitude  X,  longitude  L  and  height  h  above  the  mean  sea  level.  These  can  be  converted 
to  the  Global,  Earth-centered  coordinate  system  as  follows; 

X  =  (Re  +  h)  cos  X,  cos  L 
Y  =  (Re  +  h)  cos  X,  sin  L 
Z  =  (Re  +  h)  sin  X, 
where  Re  is  the  radius  of  the  earth. 

As  indicated  earlier,  the  North-referenced  coordinate  system  is  a  convenient 
representetion  for  the  bistatic  geometry  and  has  been  employed  in  this  work.  Next, 
we  focus  on  the  conversion  of  data  presented  in  the  Global,  Earth-centered  coordinate 
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z 


Fig.  3-5:  Global,  Earth-centered  Coordinate  System 


3 


system  to  the  North-referenced  coordinate  system. 

The  measxirements  of  the  total  range  and  baseline  are  readily  available,  and  are 
independent  of  the  coordinate  system  in  use.  Therefore,  their  conversion  need  not  be 
considered  here.  On  the  other  hand,  the  angles  ftr  and  Or  in  the  bistatic  plane  are  not 
available  as  direct  radar  measurements,  when  Global,  Earth-centered  coordinate 
system  is  used.  The  calculations  of  &r  and  0r,  using  the  data  given  in  the  Global, 
Earth-centered  coordinate  system,  are  straightforward.  Let  us  consider  the  system 

shown  in  Fig.  3-6.  Using  the  position  vectors  of  the  transmitter  =  (Xp  Yp  Zp)^,  the 

receiver  =  (Xr,  Yr,  Zr)’’,  and  the  target  =  (X^,  Y,^,  the  following 

vectors  can  be  formed: 


trr  =  h-K  ’  (Xt-X,,. 

(3.25) 

~  ^R~^Tgt» 

(3.26) 

and 

Lr  =  ?p-?R  =  (Xp-Xr,  Yp-Yr,  Zp-Zr)^ 

(3.27) 

Then,  from  the  definition  of  the  inner  product,  the  angle  between  Lpp  and  Lr  is  given 

by 

«  =  cos  ' 

■  1^1 

(3.28) 

and  the  angle  between  Lrp  and  (-Lg)  is  given  by 
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Tgt 


Fig.  3-6  :  Geometry  for  the  conversion  of  the  Global,  earth-centered  coordinate  system 
to  the  North-referenced  coordinate  system 
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(3.29) 


_  -1  ^ 

Y  =  COS  - — 

mcrl  ’  1^1 

where  A  •  fi  is  the  inner  product  of  vectors  A  and  fi  and  |Al  is  the  magnitude  of  the 
vector  A.  Since,  0  i  cos  'C)  <  the  angles  6.^  and  e^,  in  the  bistatic  plane,  can  be 
obtained  as 

— -o,  when  target  is  in  the  northern  hemisphere 

9,  =  (  ^  (3.30) 

—  +a,  when  target  is  in  the  southern  hemisphere 


— -Y,  when  target  is  in  the  northern  hemisphere 

e,  =  {  ^ 

—  +Y  when  target  is  in  the  southern  hemisphere 
2 

Using  the  above  procedure,  measiurements  given  in  the  Global  earth-centered 


coordinate  system  can  be  converted  to  the  North-referenced  system. 


3.3  T.ncal  Vertical  CoorHinafp  System 

Another  typical  situation  encountered  in  practice  is  that  the  measurements  are 
taken  in  a  local  vertical  coordinate  system  centered  at  the  receiver  site.  This  local 
coordinate  system  is  usually  specified  with  the  receiver's  local  vertical  as  the  z-axis, 
and  the  local  horizontal  as  the  x-y  plane,  in  which  the  x-axis  points  in  a  convenient 
direction,  such  as  towards  the  North  pole  and  the  y-axis  points  towards  east  This 


configuration  is  known  as  the  north-east-down  coordinate  system.  In  this  setup,  the 
target  angle  of  arrival  is  measured  in  terms  of  the  azimuth  angle  Ar,  positive 
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clockwise  from  the  x-axis^  and  the  elevation  angle  Ek,  measured  from  the  receiver's 
nadir.  Note  that  the  measurements  of  the  total  range  and  the  baseline  are 
independent  of  the  coordinate  system  in  use.  Therefore,  their  conversion  need  not  be 
considered.  In  the  following,  we  consider  the  derivation  of  6^^.  The  angle  6^  can  be 

obtained  using  the  formula  given  in  Section  3.1.  Assume  that  the  receiver  has  an 
estimate  of  the  transmitter's  angular  position  A^t  and  Ert.  Also  assume  that  the  target 
and  transmitter  are  both  above  the  local  horizontal,  Le.,  the  x-y  plane.  Other 
situations  will  be  considered  later.  Consider  the  geometry  shown  in  Figure  3-7  (case 
I).  Note  that  Tgt,  K  and  M  form  a  plane  parallel  to  the  local  horizontal  where  K  and 
M  denote  the  intersections  of  this  plane  with  the  z-axis  and  the  baseline  respectively. 
Using  the  cosine  law  in  the  triangle  Tgt  K  M,  we  have 

TgTM  2  =  2  _  2  T^tK  m  co8(A^-A^t) 

Next  consider  the  bistatic  plane,  i.e.,  the  plane  formed  by  T^  Tgt  and  R,.  Using  the 
cosine  law  for  the  triangle  with  vertices  M,  Tgt  and  R,  in  the  bistatic  plane,  we  have 


TgtM  2  =  TgtR,  2  +  R^  2  -  2  TgtR,  R^  cos(9O°+0b) 

=  T^t^2  +1^^  +2T^^sin0ji 


But 


TgtK  =  TgtRj  cos(Ej-9(3^  =  TgtR,  sinE^ 


(3.33) 


(3.34) 


R,K  =  TgtR^  sm(E|i-90®)  =  R^M  sin(ERT-90^ 


and 
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Fig.  3-7:  Geometry  for  the  conversion  of  the  local  vertical  coordinate  system  to 
the  North-referenced  coordinate  system  for  case  I 
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(3.36) 


Equating  the  right-hand  sides  of  (3.39)  and  (3.40),  and  then  dividing  by  TgtRj^  we 
have 
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(3.41) 


Rearranging  the  two  sides  of  (3.41),  we  arrive  at 

-  sin0n  =  cosEj^cosEg^] 

Note  that,  when  >  0,  the  target  is  in  the  northern  hemisphere;  and  when 

Aj-A^  <  0,  the  target  is  in  the  southern  hemisphere.  Thus, 


-sin‘^[cos(AR-A|ix)sln  Eg  sin  Ert  +  cos  Eg  cos  Erj],  when  Ar-Arx  >  0  (3.43a) 

7C  +  sin'^[cos(AR-ARx)sin  Er  sin  Erj  +  cos  Er  cos  ErxI,  when  Ar-Arx  <  0  (3.43b) 


Next  consider  the  case  where,  as  shown  in  Fig.  3-8,  both  the  target  and  the 
transmitter  are  below  the  local  horizontal,  Le.,  below  the  x-y  plane  (case  II). 
Compared  with  the  prior  situation,  it  is  obvious  that  the  derivation  of  6r  can  be 
carried  out  using  the  same  procedure.  In  fact,  the  only  modifications  needed  are  the 
substitutions  of  ji-Er  and  ti-Erx  for  Er  and  Erx,  respectively,  in  equations  (3.34),  (3.35) 
and  (3.36).  Thus,  the  resxilt  of  (3.42)  still  holds  true  in  this  case.  Finally,  consider  the 
case  when  the  target  and  the  transmitter  are  on  different  sides  of  the  local  horizontal 
of  the  receiver  (case  III).  Without  loss  of  generality,  we  assume  that  the  target  is 
above  the  plane  and  the  transmitter  is  below  it,  as  shown  in  Fig.  3-9.  Note  that  point 
T  is  the  mirror  image  of  the  transmitter  T^  with  respect  to  the  x-y  plane.  Also,  Tgt, 
K  and  M  form  a  plane  parallel  to  the  x-y  plane.  The  point  S  is  the  intersection  of  the 
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Note:  ^  T,  is  out  of  the  x-z  plane,  towards  the  viewer. 

Tgt  is  out  of  the  x-z  plane,  away  from  the  viewer. 

N  IS  the  north  of  the  North-referenced  coordinate  system 


Fig.  3-8  :  Geometry  for  the  conversion  of  the  local  vertical  coordinate  system  to  the 
North-referenced  coordinate  system  for  case  II 
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Fig.  3-9  :  Geometry  for  the  conversion  of  the  local  vertical  coordinate  system  to 
the  North-referenced  coordinate  system  for  case  III 


^3 


line  and  the  line  normal  to  the  x-y  plane  and  passing  through  M.  Note  that 

TgtM  ±  MS,  and  R^M  =  R^S.  Also,  in  the  triangle  formed  by  M,  S  and  R,,  it  is  found 
that 


MS  2  =  R^M  2  +  R^S  2  -  2  R^M  R^S  cos[2(90°-En^)]  =  2  RjjM  ^(1 -cos 2(90° 

=  2  r;^  2(1-(cos2(90°-Ejit,)  -  suiW-Ert)))  =  2  RJa  ^(l-sin^T  +  cos^E^^) 
=  4  RjjM  ^  cos^Ejj^ 

(3.44) 


where  we  have  used  the  fact  that  the  x-y  plane  bisects  the  angle  ZTR^T,,  and  the 
lower  half  of  this  angle  is  equal  to  90‘’-Erx.  However,  in  the  triangle  formed  by  M,  S 
and  Tgt,  we  have 


TgtS  2  =  TgtM  2  +  MS^ 

=  Tpi  2  +  4  2  cos^  Ejix 

where  (3.44)  has  been  used.  In  the  triangle  formed  by  Tgt,  S  and  R,,  using  the  cosine 
law,  we  have 

TgtS  2  =  T^  2  ^  R^S  2  -  2  r;;S  cos(9O°+0b) 

(3.46) 

=  Tgt^  2  +  ^;m  2  +  2  Tgts;;  r;;m  smOj^ 

where  R^S  =  R^^M  has  been  invoked.  Combining  (3.45)  and  (3.46),  and  rearranging 
terms,  we  have 

TgtM  2  =  TgtR;;  2  +  r;;m  ^  -  4  ^  ^  cos^Ej^^  +  2  Tgt^  ^  sm0ji  (3-47) 
Also,  in  the  triangle  formed  by  Tgt,  K  and  M,  using  the  law  of  cosines. 
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(3.48) 


TgtM  2  =  TgtK  2  +  KM  ^  -  2  TgtK  KM  cosCA^-A^t) 
From  the  triangle  formed  by  Tgt,  K  and  we  have 


TgtK  =  TgtR,  cos(Ej-90^  =  TgtR,  sm% 


(3.49) 


R^K  =  TgtR,  sm(E|t-90P)  =  R,M  anfOOP-E^T) 


(3.50) 


KM  =  R,M  cos(90P-E|„.)  = 


(3.51) 


where,  in  (3.50)  and  (3.51),  we  have  again  used  the  fact  that  the  x-y  plane  bisects  the 


angle  ZTR,T,.  From  (3.50), 


R,M  =  TgtR,  sin(Ea-90°)  [sm(90°-EaT)]’'  =  TgtR,  I 


Substituting  (3.52)  into  (3.51),  we  have 


(3.52) 


_  _ (  cosE. 

KM  =  TgtR. - P-  smE|tj 

cosEj^rj 

Substituting  (3.49)  and  (3.53)  into  (3.48),  we  obtain 


(3.53) 


2  f 


TgtM  ^  =  TgtR,  ^  sm%  +  TgtR^  ^  ^ 

_  ,  COSEg  »  A  X 

+  2  TgtR,  *  sinEj  - —  smE^T-cosCAj^-AnT) 

I  cosEn^l 


(3.54) 


Substituting  (3.52)  into  (3.47),  we  have 
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TgtM  ^  =  TgtR^  2  +  TgtR,  ^ 
/cosEj) 


f 


\2 


-TgtR, 


[cosBt, 


i«»V 

sme^ 


-4TgtR, 


cosE^. 


\2 


cos^  Ejj 


t; 


(3.55) 


Equating  the  right  hand  sides  of  (3.54)  and  (3.55)/  and  dividing  the  result  by  TgtR^ 
we  have 


sin^Ej 


'  cosE|^ ' 

,  (  cosE. 

sm^Ep-  +  2  sinEp - — 

cosEg^^^ 

IcosExtJ 

=  1  + 


cosE|i 


^2 


[cosEj. 


t; 


4  cos^E|^  - 


cosE^N 


smE„  co8(A^-A„) 
sm6|^ 


(3.56) 


After  rearranging  terms  in  (3.56)/  we  again  arrive  at  the  same  result  as  (3.42).  AlsO/  the 
analysis  for  the  case,  in  which  the  target  is  below  the  local  horizontal  of  the  receiver/ 
while  the  transmitter  is  above/  it  takes  exactly  the  same  steps  as  before/  and  results 
in  the  same  equation  as  (3.42).  Therefore/  it  is  concluded  that  (3.43)  properly  describes 
the  general  situation.  Then/  with  the  expression  (3.43)  and  using  the  procedures 
developed  in  Section  3.1/  9r  can  be  obtained. 

3.4  Disctission 

In  this  chapter/  we  have  desaibed  the  two-dimensional  North-referenced 
coordinate  system.  It  adequately  describes  the  bistatic  geometry.  Due  to  its 
convenience/  we  have  used  it  in  our  analysis.  We  also  described  two  other  widely 
used  coordinate  systems  and  provided  transformation  procedures  for  conversion  to  the 
north-referenced  coordinate  system.  It  should  be  noted  that  conversion  from  the 
north-referenced  coordinate  system  to  other  higher  dimensional  systems  is  not  possible. 
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4  ANALYSIS  OF  AMBIGUITY  FUNCTION  FOR  BISTATIC  RADAR 
This  chapter  investigates  the  concept  of  ambiguity  function  for  a  bistatic  radar. 
First,  we  consider  the  return  signal  from  a  moving  tai^et  for  a  bistatic  radar  and 
derive  expressions  for  target  position  and  its  velocity.  This  problem  has  been  studied 
extensively  for  a  monostatic  radar  [2]  but  no  results  for  the  bistatic  case  are  available 
and  are  developed  here.  Next,  bistatic  geometry  factors  are  included  in  the  evaluation 
of  the  ambiguity  hmction.  It  is  shown  that  plotting  the  ambiguity  function  on  the 
delay-Doppler  plane  is  no  longer  appropriate.  Suitable  arguments  for  the  ambiguity 
function  are  determined.  Examples  are  given  to  illustrate  the  concepts  developed  in 


this  chapter. 

4.1  Measurements  on  A  Moving  Target 

Consider  the  bistatic  geometry  consisting  of  the  transmitter  T^  the  receiver  R, 
and  the  target  Tgt  as  shown  in  Fig.  4-1.  Assume  that  a  signal  is  transmitted  by  T^,  at 
time  t=0.  It  reaches  the  target  at  t=Xi  and  the  reflected  signal  is  received  by  R,  at  t=x,. 
Then  the  total  time  delay  x(t),  i.e.,  the  total  travel  time  from  T,  to  R^  is  given  by 


The  time  delay  x(t)  is  a  hmction  of  the  bistatic  geometry,  i.e.,  the  relative  positions  of 
T^  Tgt  and  R^  The  same  value  of  x(t)  is  attained  for  targets  lying  on  the  same 
constant  range  contours.  Therefore,  we  may  express 


c  c 

where  c  is  the  wave  propagation  speed,  R(*)  is  the  total  range,  Rx(*)  aud  Rr(*)  are  the 
range  components  corresponding  to  the  transmitter  and  receiver  respectively.  Note 
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Tgt  ( signal  reaches  target 


Tx  ( signal  transmitted 
at  t  =  0 ) 


Rx 

( signal  received 

at  t  =  Ta ) 


TjjTgt  =  ( speed  of  light  )xxi 

TgtR,  =  ( speed  of  light  )x(  x^-Xi ) 

total  time  delay  ( total  travel  time  of  the  signal )  = 


Fig.  4-1:  Timing  relationships  in  bistatic  radar  geometry 
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Similarly  using  (4.2)/ 

x(xj  =  -R(x,-Th(t^) 
c 


(4.6) 


Using  Taylor  series  expansion  to  expand  x(t)  around  x,  we  can  write 


Let  us  next  consider  the  first  two  derivatives  of  x(t)  in  the  above  expansion.  From 
(4.6), 
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=  -R(t-Vt))[l-TR(t)]  1,^^ 

c  * 


-,(•))]  I..,. 


>  -[R(t-r,(t))Il-*»(t)P-R(t-t,(t))T,(t)]  I,., 

c  * 


where  (4.8)  has  been  used  in  deriving  (4.9).  From  (4.5) 


c  * 


(4.10) 


Rearranging  (4.10),  we  have 


Also,  using  (4.10)  and  (4.11) 


(4.11) 


Thus,  rearranging  (4.12) 
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(4.13) 


Using  (4.11)  and  (4.13),  we  can  rewrite  (4.8)  and  (4.9)  as 

lx(T^  =  ^R(t-x^(t))  [l+l6^(t-x^(t))]-M,.,^ 

and 

—  x(x^  =  M(t-xa(t))  [l+-ft»(t-x^(t))]-2-lR(t-Xa(t))  [-R|i(t-X|i(t))]  [l+-^R|i(t-x^(t)) 

<Jt2*^C  C  C  C  c 

=  -[R(t-X^(t))(l^l^(t-Xj(t))) 

c  c 

-A(t-ta(t))(lR^(t-Xj(t)))]  [1  +-^(t-x^(t))]-’  1,.,^ 

c  c 

=  -[(RT(t-'Ca(t))+R»(t-Xj(t)))(l+-^(t-Xj(t))) 

c  c 

-  (Rr(t-t,(t))+R»(t-Xa(t)))(^R^(t-x^(t)))]  [1+-R»(t-Xa(t))]-»|,.,^ 

c  c 

=  ^l^t-x^(t))+R^(t-x^(t))+lfer(t-x^(t))^(t-x^(t))+lR.(t-x^(t))^(t-Xa(t)) 

c  c  c 

-  -^(t-‘rR(t))^t-tR(t))--R»(t-T^(t))R»(t-x^(t)M  [l+X(t-CR(t))]-'lt., 

c  c  c  * 

=  -[(l+-R»(t-x^(t)))RT(t-x^(t)) 
c  c 

c  c 

(4.15) 

where,  in  deriving  (4.15),  we  have  used  the  relationships 
R(t)  =  RT<t)+R^(t)  and  R(t)  =  Rj(t)+R^(t). 

As  a  check,  if  we  set  Rj,(t)  =  Rn(t)  =  |R(t)  and  Xjj(t)  =  ^^(1),  equations  (4.6), 
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(4.14)  and  (4.15)  reduce  to  the  well  known  results  for  the  monostatic  case  [2]. 


Specifically,  (4.6)  can  be  written  as 

=  -R(T,-Ta(T^) 
c 

C  Z 

= 

c  2 

where  (4.1)  has  been  used,  i.e.,  x(t,)  =  x,.  Also,  (4.14)  can  be  written  as 

at  *  c  c 

C  i  C  Z 

2  •  "c.  1  .  X.  , 

c  2  c  2 

and  (4.15)  can  be  written  as 


(4.16) 


(4.17) 
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1  [(l*iRi(x.-t,CV))  ft/t.-t,(T^)*(l-ift,(T.-T,CT.)))  fi,(T,-T»(T^)] 
c  c  c 


c 

1  [(1  ^(1  -7ftR(^.-^'(\)»  j)i 


c  c 


c  z 


(4.18) 


For  targets  of  practical  interest,  ^(t)/c  and  R^CO/c  are  very  small.  Consequently, 


(4.14)  and  (4.15)  reduce  to 


(4.19) 


-  M(t-Tii(t))| 

dt^  c 


(4.20) 


Substituting  the  expressions  for  the  derivatives  in  the  series  expansion  (4.7),  we  have 


T(t)  =  x(tJ  + 


v2 

'C,+  (t-xJ[M(T.-TR(T.))]  +— ;^[-^R(T,-T,i(T^)]  +• 

c  ^  c 


(4.21) 


In  the  analysis  to  follow,  we  will  assmne  that  the  total  range  is  a  relatively 
smooth  fimction  of  time,  over  the  coherent  processing  period,  such  that  the  second 
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and  higher  order  derivatives  can  be  dropped,  Le.,  •••  are 

negligible.  Note  that  this  assumption  may  not  hold  true  for  a  maneuvering  target 
It  will  further  be  assumed  that  the  transmitted  radar  signal  is  narrowband,  so  that  the 
complex  envelope  of  the  signal  is  insensitive  to  the  target's  motion.  That  is, 

•  f(t-T^  W'22) 

c 

where  f(t)  is  the  complex  envelope  of  the  transmitted  signal.  In  short,  we  have 

asstuned  that  the  range  rates  of  the  second  and  higher  orders  are  negligible,  and  the 
effect  of  the  first  order  range  rate,  upon  the  complex  envelope,  can  be  ignored.  The 
validity  of  these  assumptions  has  been  generally  accepted  for  monostatic  radars  [2]. 
Yet,  little  consideration  has  been  given  to  the  bistatic  case.  From  the  analysis  in  the 
following  subsection,  it  is  found  that,  for  the  same  target  motion  relative  to  the 
ground,  the  resulting  range  rates  are  smaller  in  the  bistatic  radar  than  in  the 
monostatic  case.  Therefore,  in  this  work  these  assumptions  are  employed  for  the 
ambiguity  function  analysis  of  bistatic  radar.  However,  a  detailed  study  of  this  issue, 
Le.,  the  conditions  under  which  these  asstunptions  are  valid,  is  desirable  and  may  be 
carried  out  in  the  next  phase  of  this  effort 

With  the  aforementioned  simplifications,  the  retiun  signal  can  be  modeled  as 

f(t)  =  6^(t-T(t)) 

=  6^t-T^  OstsT  (4.23) 

where  6  is  given  by  the  reflectivity  of  the  target,  E,  is  the  energy  of  the  transmitted 
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signal,  cOo  is  the  carrier's  angular  frequency,  T  is  the  duration  of  the  signal,  andfift) 

is  the  complex  envelope  of  the  additive  noise. 

Note  that  the  term  x,  in  (4.23)  is  given  by  the  total  range,  and  the  term 

represents  the  Doppler  shift  Then,  the  delay  x,  and  the  Doppler  shift 

can  be  extracted  by  using  a  matched  filter  and  a  square-law  envelope 
detector,  as  described  in  Chapter  2. 

In  the  following,  we  will  examine  the  manner  in  which  the  measurements  of 
range  and  range  rate  are  related  to  target  location  and  its  velocity.  That  is,  for  a 
hypothesized  target  location  and  its  velocity,  the  resulting  total  time  delay  and 
Doppler  shift  will  be  derived. 

4.1.1  Range  and  Target  Location 

In  bistatic  radars,  the  target  location  may  be  specified  either  with  respect  to  the 
receiver,  or  the  transmitter.  That  is,  in  the  North-referenced  coordinate  system,  the 
range  Rr  and  the  angle  9r,  or  the  range  Rj  and  the  angle  67,  are  given  to  specify  the 
target  location.  Given  the  baseline,  the  calcxilation  of  the  resulting  total  range  from 
a  given  target  location  is  then  straightforward. 

In  the  case  where  target  location  is  given  in  terms  of  Rr  and  0r,  recall  equation 

(3.1), 

R  =  Rj+^R^+L^+2R^Lsm^ 

Thus,  with  the  additional  knowledge  of  the  baseline,  the  total  range  is  determined. 
To  emphasize  the  functional  relationship  between  the  total  range  and  the  parameters 
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Rji,  L  and  6k/  we  define  the  function 


Now,  the  total  time  delay,  or  equivalently  the  total  range,  can  be  found  using  the 
hmction  D.  In  other  words,  when  bistatic  geometry  is  given  in  terms  of  Rr,  L  and  6r, 
or  in  terms  of  Rj,  L  and  0x,  the  resxilting  total  time  delay  is  given  by  Dg,  or  Dj, 
respectively. 


Ranee  Rate  and  Target  Velocit 


As  derived  earlier,  the  Doppler  shift  is  given  by 


Since  the  time  argument  in  the  measurements  is  not  of  concern  in  this  analysis,  it  will 
not  be  expressed  explicitly  in  the  following  derivation.  We  will  first  consider  the 
generalized  case  where  R,  and  Tgt  are  all  moving,  as  shown  in  Fig.  4-2,  and  then 
focus  on  the  effects  of  target  motion  alone.  Note  that  the  angle  (|)t  is  the  angle 
between  the  transmitter's  velocity  vector  Wj  and  the  north  of  the  North-referenced 
coordinate  system,  while  the  angle  (t>R  is  the  angle  between  the  receiver's  velocity 
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vector  Vr  and  the  north.  Both  ())x  and  (t)R  are  measured  in  the  same  convention  as  6t 
and  0R,  i.e.,  they  are  positive  clockwise.  The  angle  (])  is  the  angle  between  the  target 
velocity  vector  V  and  the  bistatic  bisector,  positive  clockwise  from  the  bisector.  By 
definition,  it  can  be  written  that 


dt^  dt^ 


From  Fig.  4-2,  it  is  easy  to  see  that 


(4.27) 


i-Rj.  =  -V  cos  cos(<t)x-eT) 

dt  2 

Note  that  the  first  term  on  the  right  hand  side  of  (4.28)  is  the  velocity  projection  of 
V  along  the  line  connecting  Tgt  and  Tx,  and  the  second  term  the  projection  of  Vj 
along  the  same  line.  Also,  we  find  that 


— Rj^  =  -V  cos(<j)+-L&J-)-Vji  cos(<j>R-0ii)  (4.29) 

dt  2 

Here,  in  (4.29),  the  first  term  on  the  right  hand  side  represents  the  velocity  component 
of  V  along  the  line  connecting  Tgt  and  Rr,  and  the  second  term  the  component  of  Vr 
along  the  same  line.  The  appearances  of  the  absolute  sign  in  (4.28)  and  (4.29)  are  due 
to  the  way  p  is  defined.  Recall  that  P  =  ej-QR.  Then,  0  <  p  ct  when  the  target  is  in 
the  northern  hemisphere,  and  -jt  <  P  <  0  when  the  target  is  in  the  southern 
hemisphere. 

Thus,  substituting  (4.28)  and  (4.29)  into  (4.27),  we  have 
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— R  =  -V[C08((|)-^)+  C08(4>+-^)]-Vt  C08(^-eT)-V^  COS(<|)R-e|l) 
dt  2  2 


=  -V[cos4>cos— +sin4)siii-^^ +  cos4)COs— -sin4>sin-^^] 
2  2  2  2 


-  V^cosC^t-Bt)  -  VnCOs(^j-eg) 


=  -2Vcos4)COs-^ -VtCOs(^-0t)-VhCO8(^-0ii) 


(4.30) 


The  last  two  terms  on  the  right  hand  side  of  (4.30)  are  due  to  transmitter  and  receiver 
motion.  In  this  project,  we  assumed  that  these  two  quantities  are  known,  Le., 
transmitter  and  receiver  are  operating  in  a  cooperative  mode  and  the  Doppler  due  to 
their  motion  can  be  subtracted  from  the  measured  Doppler.  Equivalently,  it  was 
assumed  that  both  the  transmitter  and  the  receiver  are  stationary  in  the  following 
analysis.  Now,  the  remaining  term  is  the  one  due  to  target  motion.  This  term 
indicates  that  the  resulting  Doppler  shift  is  determined  by  the  projection  of  the  target 
velocity  along  the  bistatic  bisector.  Note  that  the  Doppler  shift  depends  on  the  target 

location  in  the  bistatic  geometry  and  it  appears  in  (4.30)  via  the  term  cos-^.  Also,  due 


to  the  mriltiplication  factor,  cos  1,  the  Doppler  shift  due  to  the  target  motion  for  the 

bistatic  radar  is  always  smaller  than  or  equal  to  |2V  cos  (1)|,  the  corresponding 
Doppler  shift  for  the  monostatic  case. 

Next,  we  derive  an  expression  for  cos-|  in  terms  of  the  bistatic  geometry 
parameters  L,  Rr,  and  9r.  From  trigonometry 


59 


_  6  1+cos  B 

cos  *  - - 

2  ^  2 


(4.31) 


From  the  law  of  cosines. 


=  R?  +  Ri-2IlrR*  cos  p 


Rearranging  the  two  sides  of  (4.32),  we  have 


(4.32) 


cos  P  = 


2R^ 


(4.33) 


Substituting  Rj.  =  jRi+L^+2Rj^Lsm0n  into  (4.33),  we  have 


cos  p  = 


2R^-t-2R|^Lsin6|^ 


Ri(+LsmB|^ 


2R^jR^+L^+2R^Lsin0a  Lsin0j 


(4.34) 


Then,  substituting  (4.34)  into  (4.31),  we  have 


B  1  Rb+LsuiOb 

2V  cos4>cos—  =  2Vcos^  — + —  ^ 

^  2t/Ri+L2+2RgLsme^ 


(4.35) 


Again,  to  emphasize  the  functional  relationship  between  the  resulting  range  rate  and 
the  location  parameters  for  the  bistatic  geometry,  we  define  the  function 


=  W(Vcos4>,RR.L,0R)4-!i  [  2  Vcos*  -  + 


I  Ri^-^Lsin6|^ 


(4.36) 


2M+L^+2RoLsM 


That  is,  the  Doppler  shift  is  expressed  as  a  function  of  V  cos  (|>,  Rr,  L  and  Br. 


Alternatively,  substituting  R^  =  ^Rj+L^-2RjLsm0.j.  into  (4.33),  we  have 
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COSp  = 


2Il?-2RjLsineT 


Rj.-Lsin0T 


(4.37) 


211^/11? +L^-2RrLsmeT. 

Then,  substituting  (4.37)  into  (4.31),  we  have 

2Vcos*cosi  =  2VCOS*  (4.38) 

^  \j  ^  2^J^-2RjLsmeT 

We  may  define  the  Doppler  shift  Wt  as  a  function  of  V  cos  Rj,  L  and  0t,  as  follows 


W.  =  W(Vcos<|),Rt,L,-0t) 


(i).  1  R#m  Lsiii0>«i 

=  -^[2Vcos<|)  -  +  ^  — 


(4.39) 


2  Jr?  +  L2-2Ri.Lsin0T 


Thus,  with  (4.36)  or  (4.39),  the  resulting  Doppler  shift  is  described  in  terms  of  the 


target  location  for  the  bistatic  geometry. 

In  the  analysis  thus  far,  we  have  derived  the  target  retiun  model  for  bistatic 
radars.  The  geometry  factors  in  the  measurements  have  also  been  examined  by 
relating  the  measurements  with  locations.  In  the  following,  we  will  consider  the 
consequences  on  the  representation  of  ambiguity  function  for  bistabc  radars. 


4.2  Ambieuitv  Function  Plots  for  Bistatic  Configuration 


In  radar  systems,  waveforms  are  selected  based  upon  their  ability  to  satisfy  the 
requirements  for  detection,  measvirement  accuracy,  resolution,  ambiguity  and  clutter 
rejection.  Ambiguity  function  plots  are  examined  for  a  qualitative  determination  of 
the  suitability  of  different  waveforms  in  meeting  the  above  system  requirements  [1]. 
In  practice,  ambiguity  function  is  plotted  as  a  function  of  delay  and  Doppler  derived 
from  the  return  signal.  In  monostatic  radars,  the  use  of  delay  and  Doppler  as  the 


arguments  of  the  ambiguity  function  are  adequate  because  of  their  linear  relationship 
with  the  desired  quantities,  Le.,  target  range  and  the  range  rate.  In  particular,  the 
delay  is  equal  to  two  times  the  range  divided  by  the  speed  of  light  and  Doppler  is 
equal  to  two  times  the  range  rate  divided  by  the  wavelength.  Therefore,  the  pair 
delay  and  Doppler  or  the  pair  range  and  range  rate  may  be  used  in  an  interchangeable 
manner  as  the  arguments  of  the  ambiguity  function.  Note  that  once  the  target  range 
and  range  rate  are  known,  its  position  and  its  relative  (radial)  velocity  with  respect 
to  the  monostatic  radar  are  known  completely. 

For  a  bistatic  radar,  the  above  observations  are  no  longer  valid  due  to  the 
geometry  of  the  bistatic  configuration.  The  transmitter  and  the  receiver  are  not  at  the 
same  site  and,  therefore,  the  relative  position  and  the  relative  velocity  of  the  target 
needs  to  be  determined  with  respect  to  either  the  transmitter  or  the  receiver  or  some 
other  suitable  reference  point  In  addition,  the  effect  of  the  transmitter  and  the 
receiver  not  being  at  the  same  site  should  be  taken  into  account  while  computing  the 
ambiguity  function.  In  Section  4.1,  two  fimctions  Dt(')  and  Wj(*)  with  respect  to  the 
transmitter  and  fimctions  Dr(0  and  Wr(*)  with  respect  to  the  receiver  were  defined. 
Here,  we  focus  our  attention  on  the  evaluation  of  the  ambiguity  function  with  respect 
to  the  receiver.  Similar  results  can  be  obtained  with  respect  to  the  transmitter  or  any 
other  reference  point 

In  the  bistatic  configuration,  the  baseline  L  is  assumed  to  be  given.  When  we 
evaluate  the  ambiguity  function  with  respect  to  the  receiver,  the  receiver's  look  angle 
Bn  can  also  be  assumed  known.  The  actual  position  of  the  target  can  be  determined 
from  the  knowledge  about  the  distance  between  the  target  and  the  receiver,  Le.,  the 
range  Rr..  The  total  time  delay  based  on  these  system  parameters  and  the  bistatic 
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geometry  can  be  expressed  as 


Dr,  =  -[Rr,  +  2R^  L  sin 

Also,  the  Doppler  shift  due  to  target  motion  which  takes  into  account  the  bistatic 
geometry  can  be  expressed  as 


0>r 


W„  =  —  [2V.  cos4> 
^  c 


1  . 


+  L  sin  0^ 


2^  +  +  2R^  L  sin  Or 


where  V,  is  the  target's  actual  velocity  along  the  bistatic  bisector.  As  in  the 
monostatic  case,  the  ambiguity  function  for  the  bistatic  radar  can  be  expressed  as  a 

function  of  the  delay  and  Doppler,  Le.,  as  6CDr^4Jr,.Wj^,Wj^).  As  discussed  before, 

the  delay  and  Doppler  terms  include  the  effects  of  bistatic  geometry.  Note  that  the 
subscript  "H"  denotes  the  delay  and  Doppler  bin  that  is  under  investigation.  The 
peak  of  the  ambiguity  function  occtus  at  the  point  where  Dj^=Dj^  and 

While  the  ambiguity  function  expressed  in  terms  of  Dr  and  Wr  provides  some 
information,  its  utility  is  limited.  This  is  because  in  most  applications,  resolution  in 
terms  of  the  position  and  velocity  of  the  target  is  desired.  Due  to  the  nonlinear 
relationship  between  the  pair  delay  Dr  and  Doppler  Wr  and  the  pair  Rr  and  V,  one 
cannot  immediately  observe  the  behavior  of  the  waveforms  as  a  function  of  Rr  and 
V.  This  is  especially  inconvenient  due  to  the  fact  that  the  primary  purpose  of  the 
ambiguity  function  is  to  provide  a  visual  representation  of  waveform  behavior  as  a 
function  of  position  and  velocity.  Therefore,  we  suggest  the  use  of  Rr  and  V  as  the 
arguments  of  the  ambiguity  function.  This  wiU  enable  one  to  examine  the  suitabiUty 
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of  different  waveforms  in  a  convenient  manner.  In  the  following,  we  present  two 
examples  where  we  consider  a  Gaussian  pulse  and  a  pulse  train  consisting  of  three 
rectangular  pulses.  Ambiguity  function  plots  as  a  function  of  Rr  and  V  for  a  variety 
of  bistatic  configurations  are  given.  These  exhibit  the  effects  of  geometry  on  the 
ambiguity  function.  In  these  figures,  Rr  represents  Rr,  and  both  R„  and  Rr,  represent 

Example  1:  Simple  Gaussian  Pulse 

We  first  consider  a  transmitted  waveform  whose  complex  envelope  is  a  simple 
Gaussian  pulse.  Let 


f(t)  =  (-^)^  e 


■ooit;£“> 


(4.40) 


Note  that  f(t)  has  unit  energy  since 


|f(t)P  =  (^)’  e 


(4.41) 


IS  a  zero  mean  Gaussian  function  with  variance  y .  The  target  return  can  be  written 
as 


f(t-t )  e'"'-*''-’  =  (-5-)i  e 

The  time-frequency  autocorrelation  function,  as  defined  in  Chapter  2,  is  given  by 
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4>('Ch»^4»“Dh'“d,)  =  /  f  (t  -Ch)  ®  * 


(4.43) 


Let  z=t-^.  Then, 
2 


.  V  /•  +  fVr-XZl)  gJ“o/**"T"^  ■  2 

♦(Tip^.>“Dh'"d.)  =  J  f(z+— )  f  (z  ^  )  e 


^  J  f(z+:!»fi)  f*(z-^)  e^‘^‘"”‘^dz  (4-44) 


Substituting  (4.40)  into  (4.44),  we  have 


^K«-i>h"“d^— >  f  /  1  e 

=  ®  J  V  ^ 


1  (.♦3:i)*/rr»  (z-^f/zr^ 
i-)*  e  *  e  * 


(4.45) 


Completing  the  square  in  the  exponent  in  the  integral,  and  simplifying,  we  have 

A.  'I  -  f  (-Ly  e’  ^ 

—00 


•  e  ^  “  dz 

_ _ 


(4.46) 


where  the  identity 


i~y  /  6  <i 


has  been  used.  Then,  by  definition,  the  ambiguity  function  is  given  by 


_  a  21*  * 


(4.47) 


Figure  4-3  is  the  plot  of  the  ambiguity  function  in  the  delay-Doppler  plane  given  by 

(4.47)  with  parameters  7=1.2x10“*  sec,  x^=4xl0"*  sec  and  =1200  latVsec..  For  the 

monostatic  case,  the  ambiguity  function  is  commonly  plotted  in  this  manner.  For  the 
bistatic  case,  if  the  bistatic  geometry  is  ignored  and  the  ambiguity  fimction  is  plotted 
in  terms  of  total  delay  tn  arid  Doppler  oIdh,  the  resulting  ambiguity  function  will  be 
the  same  as  shown  in  Fig.  4-3.  Also  given  in  Fig.  4-3  is  the  equal-height  contour  plot, 
where  the  magnitude,  as  a  fraction  of  the  peak  value,  is  shown  alongside  the  contour. 

As  discussed  before,  for  monostatic  radars,  the  ambiguity  function  given  in 

(4.47)  can  also  be  expressed  in  terms  of  the  range  and  radial  velocity  as 


®  )  =  0  - -) 

HI  f*  ^  f* 


(4.48) 


where  is  the  one-way  range  and  is  the  radial  velocity  of  the  target  For 
comparison.  Fig.  4-4  is  the  plot  of  (4.48)  by  setting  the  carrier  frequency 

Wj,  =  3x10*  rad/sec,  the  actual  distance  ~  ^  ~  ^*10^  rn,  and  the  actual  radial 
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velocity  =  600  m/sec.  Note  that,  except  for  constant  scaling  factors.  Fig.  4-3 
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T=  1.2e-4(sec) 


and  Fig.  4-4  have  the  same  shape.  This  was  expected  due  to  the  aforementioned  linear 
relationships. 

However,  as  proposed  earlier,  ambiguity  function  for  a  bistatic  radar  with 
respect  to  the  receiver  needs  to  be  written  as 

transmitted  waveform's  response  to  the  target  can  be  directly  related  to  target's 
position  and  its  velocity.  In  Figures  4-5  to  4-11,  we  give  the  plots  of 
0(D  ,  D  ,  Wo  ,  W- )  for  various  values  of  receiver  look  angles.  For  these  results,  it 

has  been  assumed  that  =  6x10^  m  and  L  =  10*  m.  Also,  the  actual  velocity 

component  along  the  bistatic  bisector,  is  set  at  600  m/sec.  With  these 
specifications,  the  target  locations  considered  represent  an  operating  region  close  to 
the  receiver.  From  the  resulting  ambiguity  diagrams,  it  is  obvious  that  the  effects  of 
geometry  factors  are  more  prominent  as  the  target  approaches  the  baseline.  Moreover, 

when  the  target  is  on  the  baseline,  the  resulting  delay  is  ^  and  the  Doppler  is  zero, 

regardless  of  the  position  on  the  baseline  and  the  velocity  of  the  target  Therefore, 
the  ambiguity  function  diagram  of  equal  height  in  the  region 
0  i  Rn^  s:  1  and  -«  <  Vh  <  »  is  obtained,  when  the  receiver  look  angle  is  equal  to 

Next,  we  consider  the  target  at  a  distance  larger  than  the  baseline,  by  changing 

2 

L  to  5x10*  m.  That  is,  the  region  being  considered  encompasses  both  the  transmitter 
and  the  receiver,  Le.,  the  cosite  region.  Figures  4-12  to  4-18  are  the  plots  of  the 
corresponding  ambiguity  fimction  6(0^^,  W^^,  for  various  values  of  receiver 
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( Va=  6{)0(ni/sec) ) 


case  considered: 

0r=9O° 

L=  100km 

RRa=60km 

Va=600m/sec 


xlO'*  T=  1.2e-4(sec) 


Fig.  4-5: 


Rr  (  Rra=  6.0e4<m) ) 

Bistatic  ambiguity  function  for  the  Gaussian  pulse  for 
considered 
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8 


Rr  (  Rra=  6.0e4{m) ) 


TnT5  fTiiTB  til#  i 


( Va=  600(m/sec) ) 


case  considered: 

L=  100km 
^Ra'^^Okm 

Va=600m/sec 


> 


2.5 

2 

1.5 

1 

0.5 
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-0.5 


-1 


-1.5 


xl04  T=  1.2e-4(sec) 


2  4  6  8  10 


Rr  (  Rra=  6.0e4(m)  ) 

Fig.  4-7:  Bistatic  ambiguity  hmetion  for  the  Gaussian  pulse  for  the  case 

considered 
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( Va=  600(ni/sec) ) 


xlO'*  T=  1.2e-4(sec) 


Rr  (  Rra=  6.0e4(m)  ) 

Fig.  4-8:  Bistatic  ambiguity  fxmction  for  the  Gaussian  pulse  for  the  case 

considered 


( Va=  6{)0(m/sec) ) 


case  considered: 


eR=-60° 


xlO**  T=  1.2e-4(sec) 


Fig.  4-9:  Bistatic  ambiguity  function  for  the  Gaussian  pulse  for  the  case 

considered 


( Va=  600(m/sec) ) 


( Va=  600(m/sec) ) 


( (035/01)009  =BA  ) 


case  considered: 


(  Va=  6(X)(rn/sec) ) 


xl04  T=  1.2e-4(sec) 


Fig.  4-13:  Bistatic  ambiguity  function  for  the  Gaussian  pulse  for  the  case 

considered 


( Va=  600(m/sec) ) 


case  considered: 

eR=o° 

L=50km 

RRa=60km 

Va=600 m/sec 


xiO'*  T=  1.2e-4(sec) 


Fig.  4-14:  Bistatic  ambiguity  fxmction  for  the  Gaussian  pulse  for  the  case 
considered 


( Va=  600(m/sec) ) 


case  considered:- 


Fig.  4-15:  Bistatic  ambiguity  function  for  the  Gaussian  pulse  for  the  case 
considered 
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(  Va=  600(m/sec) ) 


xio4  T=  1.2e-4(sec) 


Fig.  4-16:  Bistatic  ambiguity  fxmction  for  the  Gaussian  pulse  for  the  case 

considered 


8' 


case  considered: 


0r=-9O° 

L=50km 

RRa=60km 

Va=600m/sec 


look  angles.  When  this  case  is  compared  to  the  previous  one  (with  larger  baseline) 
the  degree  of  resolution  lost  due  to  geometry  factors  is  less.  For  example,  compare 
Figures  4-10  and  4-17  for  0r=-85“  where  we  observe  this  effect  quite  clearly.  This 
phenomenon  is  due  to  the  fact  that  in  the  latter  case  the  total  range  R=Rt+Rr  is  much 
larger  than  the  baseline  L.  As  R  increases  with  respect  to  L,  the  degree  of  resolution 
lost  due  to  geometry  factors  will  continue  to  diminish  and  the  bistatic  system  will 
behave  more  and  more  like  a  monostatic  system.  In  fact,  in  the  limit  as  L->0,  the 
bistatic  system  becomes  a  monostatic  system. 

Example  2;  Three  Rectangular  Pulses 

Here,  we  first  consider  a  sequence  of  (2n+l)  rectangular  pulses  having  a 
constant  repetition  rate  of  1/Tp  with  a  subptilse  duration  of  T  and  then  specialize  the 
result  to  the  case  of  three  pulses.  The  duration  of  the  pulse  train  is  T^  =  2nTp+T.  In 

many  practical  situations,  T«Tp  where  Tp  is  not  necessarily  a  multiple  of  T.  Also 
note  that  duty  cycle  T/Tp  <  50%  is  the  only  case  of  interest.  Let  the  complex  envelope 
of  a  subpulse  be  denoted  by 

u(t)  =  I  ^  2  (4.49) 

[  0  ;elsewheie 

The  complex  envelope  of  the  pulse  train  can  be  written  as 

f(t)  =  — E  u(t-kTp)  (4-50) 

y^2n+l  k=-ii 

where  f(t)  has  been  normalized  to  contain  imit  energy.  Assume  that  the  target  does 
not  fluctuate  in  the  period  Tj  during  which  it  is  illuminated  by  the  signal. 
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We  first  evaluate  <|)(i:h,  x^,  Qjj)  ^  From  (4.44),  we  have 


f  f’(t-^) 

J  0  9 


1_  2  /■“.•it'-i’.-’.!’ J(,.„  .illiil) 

1  +  1  1 _  JkT.4(T-lT„-Tj)  p  2 


2n+l  k=_n  Jw, 

•  u’(t-kT 

'  P  2 


p-ifr-lTa-tj) 


(4.51) 


Let  z  =  t-kT  .  Then, 


-'H-'.,  n 


.  u(z+.i!«:?i!)  uVz-JZsIIil)  dz 


1  J 

2n+l 


(:^)(t-|th-x.|)J 
(“d  -"d)T 


J_  ^l(o.o„^.o.)(:ip)  sm[(^)(n4)Tj  ^ 
2“*'  sill((0>B^-<.>B_)i) 


(^)(t-|t„-.J) 


(“Dh^^d.)^  (4.52) 


Now  consider  the  case  Xjj-x,  >  T.  Note  that  there  is  no  overlap  of  f(t+  and 


P(t_^^)  in  (4.51),  until  Xjj-x,  =  T  -  T.  When  i^j-x,  =  T  -  T,  the  term  f(t-^^) 


is  displaced  to  the  left  by  and  the  term  f(t+-!^S-^)  to  the  right  by  Then  one 

2  2  2 


less  pulse  is  included  in  the  overlap,  as  illustrated  in  Fig.  4-19.  We  thus  conclude,  for 
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*'"h.^.“d..“d.)  -  ^  e  J  e 

/•B-,*j<r-l<,-».-T,l)  ^  l'«~’!3i)  u-(t^kT  -iZslii!)  e'*"^'"”-*  dt 

Jh-j-IcT-Wh-^-TpI)  P  2  P  2 


1  ^  gj(“DH'“l>, 


2n+l 


,,1T,  ^ 


(^^)(T-|t„-x:.-T,|) 


k=-n+l 

/ 


(“d  -“d)T 


1  ^i(“d„*<-d;(^)  ^ 


2n+l 


-“d)T 


(4.53) 


s^(“dh-“d.)t) 

Carrying  on  in  this  manner,  there  are  (2n+l)  expressions  similar  to  the  ones  given  in 
(4.52)  and  (4.53).  Thus  the  general  expression  for  the  time-frequency  autocorrelation 
function  is 


1  ^  2sin|(^)Cr-bH-.-l^pl) 


2n+l 


^  '  E  - 

k=-n 


(“d  -Wd)T 


H  *'• 

T. 


rect 


5m((Qp,-(0„X2n^l-|k|)f 


(4.54) 


Finally,  the  ambiguity  fimction  is  given  by  the  magnitude  square  of  (4.54).  That  is. 
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1_  I  2sta((^)(T-|TH-T.-ia',|)) 

(“d 


2n+l 


E 

k=-n 


T„-T-kT„  sin((o)jj -a)jj)(2n+l-|k|)^ 

rect(-5^— ! - ?)  — i  ” — I? - L 

2T 


n  1  2siiui 

=  E  l-i - 11 


E 

k=-n  2n+l 


(^^)(T-|T:„-T.-kT  I) 


(“d  -“d)T 


I- 

-a>j3X2n+l-l 

.l)i) 

sii)| 

(  T  ' 

'(0)d 

f  W  2, 

1 

„  T„-T,-kT„ 

p  rect(  »  P)  (4.55) 


where  we  have  assumed  T  <  |  Tp.  In  Figure  4-20,  we  give  the  plot  of  (4.55)  assuming 

the  carrier  frequency  Qj,  =  3x10®  rad/sec,  the  number  of  pulses  transmitted 

N  =  2n+l  =  3,  the  subpulse  duration  T  =  4x10'*  sec,  and  the  period  Tp  =  10"*  sec.  Also, 
it  is  assumed  that  the  actual  delay  x,  =  4x10"^  sec,  and  the  actual  Doppler 

ti>n  =  1200  rad/sec. 

For  monostatic  radars,  Figme  4-21  shows  the  result  of  (4.55)  plotted  on  the 
i-6v  it  is  plotted  with  the  hypothesized  one-way  range  and  the  radial 

velocity  as  variables.  Note  that,  here,  R..  =  —  =  6x10^  m  andV.,  =  —  =  600  m/sec.. 

^.2  M.  2co„ 
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T=  .4e-4(sec);  Tp=  le-4(sec) 


xl04 


One-way  Range  (  Actual  Range  =  6.0e4(m) ) 

Fig.  4-21  :  Ambiguity  function  for  monostatic  radar  in  the  range-velocity  plane 
for  the  three  rectangular  pulses  waveform 


Again,  as  before,  the  ambiguity  function  for  the  bistatic  case  will  be  the  same  as 
shown  in  Fig.  4-20  if  bistatic  geometry  is  ignored.  On  the  other  hand,  when  this  three 
rectangular  pulse  train  is  employed  in  a  bistatic  radar,  the  resulting  ambiguity 
function  with  respect  to  the  receiver  can  be  examined  by  writing  (4.55)  in  the  form  of 
0(D„  ,  ,  W„  ,  ).  The  results  are  plotted  in  Figures  4-22  to  4-26  assuming  that 

=  6xlOV  =  600  m/sec  and  L  =  10®  m,  for  various  values  of  receiver  look 

angles.  In  order  to  compare  the  ambiguity  fxmction  for  the  above  receiver-centered 
region  with  the  cosite  region,  we  change  the  value  of  L  to  5x10^  m.  The  results  are 
plotted  in  Figures  4-27  to  4-31,  for  various  values  of  receiver  look  angles.  Now,  it  is 
again  seen  that  the  geometry  factors  have  more  profound  effects  in  the  receiver- 
centered  region  than  in  the  cosite  region.  Moreover,  within  the  same  operating 
region,  the  effects  of  geometry  factors  are  more  prominent  when  the  target  is  closer 
to  the  baseline. 

4.3  Discussion 

In  the  previous  section,  it  has  been  indicated  that  the  ambiguity  function  for 
a  bistatic  radar  needs  to  be  plotted  with  respect  to  an  appropriate  point.  For  example, 
the  reference  point  can  be  the  transmitter  site,  the  receiver  site  or  some  other  point 
of  interest.  In  the  examples,  we  focussed  on  ambiguity  function  plots  with  respect 
to  the  receiver  only.  Similar  results  can  be  obtained  in  other  cases.  Two  examples 
were  considered  which  demonstrated  the  fact  that  bistatic  geometry  plays  an 
important  role  in  ambiguity  function  analysis.  It  was  observed  that,  in  the  regions 
close  to  the  receiver,  geometry  factors  have  a  more  pronounced  effect  on  the  resolution 
capabilities  of  the  transmitted  waveform  than  in  the  cosite  region.  In  fact,  resolution 
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( Va=  6(X)(m/sec) ) 


case  considered: 


0r=9OO 
L=  100km 

RRa=60km 

\a=600m/sec 

T  ^ 


xlO^  T=  .4e-4(sec);  Tp=  le-4(sec) 


Fig.  4*22:  Bistatic  ambiguity  function  for  the  three  rectangular  pulses  waveform  for 
the  case  considered 
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case  considered: 


eR=45° 

L=  100km 
Rl^=60km 


V,i=600m/scc 


y 


( Va=  600(m/sec) ) 


23456789  10 


Rr  (  Rra=  6.0e4(m) )  xlO** 

Fig.  4-24:  Bistatic  ambi^ity  function  for  the  three  rectangular  pulses  waveform  for 

the  case  considered 
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( Va=  600(in/sec) ) 


case  considered: 

0r=9O° 

L=50km 

RRa=60km 

Va=600m/sec 


Xl04 


Rr  (  Rra=  6.0e4(m) ) 

Bistatic  ambiguity  function  for  the  three  rectangular  pulses  waveform  for 
the  case  considered 


T=  .4e-4(sec);  Tp=  le-4(sec) 


Fig.  4-29:  Bistatic  ambiguity  function  for  the  three  rectangular  pulses  waveform  for 

the  case  considered 
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( Va=  600(m/sec) ) 


( Va=  600(m/sec) ) 


case  considered: 


xiO'*  T=  .4e-4(sec);  Tp=  le-4(sec) 


Fig.  4-31:  Bistatic  ambiguity  function  for  the  three  rectangular  pulses  waveform  for 

the  case  considered 
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is  totally  lost  if  the  target  is  on  the  baseline.  Note  that  these  effects  cannot  be 
observed  if  the  ambiguity  fimction  is  plotted  using  delay  and  Doppler  as  argtunents 
in  a  conventional  manner  disregarding  the  geometry.  Note  further  that  even  for  a 
simple  three  rectangular  pulse  waveform,  the  bistatic  geometry  significantly  affects 
the  ambiguity  fimction  and  the  associated  parameter  estimates.  This  effect  will  be 
even  more  pronounced  with  more  complex  waveforms  based  on  PN  sequences.  The 
methodology  presented  here  must  be  used  to  assess  the  capabilities  of  transmitted 
waveforms  for  bistatic  radar  scenarios. 
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5  CONTINUOUS-TIME  CONTINUOUS-FREQUENCY  AMBIGUITY 
FUNCTION  OF  A  SAMPLED  SIGNAL 


In  our  work  thus  far,  we  have  assumed  that  continuous  signals  are  processed 
at  the  receiver  input  In  practice,  however,  most  modem  radars  sample  the  signal 
after  it  has  been  frequency  translated  to  baseband.  This  sampled  discrete-time  signal 
is  processed  using  digital  signal  processing  techniques.  Next  we  consider  ambiguity 
function  analysis  based  on  the  sampled  signal.  First  we  consider  the  case  where 
X  ,  I ,  CO  and  co„  are  allowed  to  assume  all  possible  values.  In  this  case,  the 

H'  a'  Djj  Da 

ambiguity  function  is  a  continuous  fvmction  of  time  and  frequency  even  though  the 
received  signal  has  been  sampled. 

5.1  Likelihood  Ratio  Test  for  a  Disaete  Signal 

We  return  to  the  detection  problem  as  discussed  earlier,  but  assume  that  the 

complex  envelope  of  the  received  pulse,  f(t),  is  represented  by  N  complex  samples 
fp,  fj,...,  fjj  y  The  detection  problem  is  now  characterized  by  the  two  hypotheses 


Hr 

H<,:  f^=n^  k=0,l,...N-l 


(5.1) 


where  s^,  Sj,...,  Sjj.j  are  samples  of  the  received  signal  complex  envelope 
s(t)=  6  at  times  t  =  +  k^,  and  nj,...,  represent  samples  of  the 

complex  envelope  of  the  Gaussian  white  noise  process  n(t).  As  before,  assume  that 
6  remains  fixed  over  the  duration  of  the  signal  pulse  and  is  statistically  independent 
of  the  noise  samples.  Let  f  be  the  (column)  vector  whose  elements  are 


103 


fy  k=0,  We  introduce  the  notation 


dio  =  E[f  IHq]  and  m,  =  E[f  iHj]  (5-2) 

for  the  conditional  mean  vectors  of  f  under  hypotheses  Hq  and  Hi,  respectively.  Also, 
we  denote  the  covariance  matrices  under  Ho  and  Hi  by 

Kj  =  E[(f-m^(f-m,)»|Hi]  i=0,l,  (^-3) 

where  the  superscript  H  denotes  the  conjugate  transpose  (Hermitian).  We  may  then 
write  the  conditional  probability  densities  of  f  under  i=0,l,  as 


exp  {-(R”-mf)K,‘'(R-ng} 


i  =  0,1.  (5.4) 


Here,  [KJ  denotes  the  determinant  of  Kj.  The  likelihood  ratio  test  is  given  by 


A(R)  A 


KqI  e 
KJ  e 


Hi 

> 

< 

Ho 


Taking  logarithms,  we  obtain 


(R''-ino”)Ko‘(R-mo)  -  (R”-mf)K-'(R-mi) 
Hi 

^  to  r|  +  fin  IKJ  -  to|Ko|  A  Y  • 
Ho 


(5.5) 
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5.1.1  Zero-mean  Case 

Let  the  mean  vectors  on  the  two  hypotheses  be  equal  to  zero.  Then 

Aq  =  Aj  =  0. 

The  likelihood  ratio  test  reduces  to 


Hi 

R""  {Ko‘'  -  R  ^  Y  . 

Ho 

Denote  the  inverse  of  the  covariance  matrix  Kj,  i=0,l,  by  Q, 

Qj  =  Kj'^  i=0,l. 

The  likelihood  ratio  test  can  then  be  written  as 


(5.6) 


(5.7) 


Hi 

R“(Qo-Qi)R  <  Y. 

Ho 

5.1.2  Zero-mean  Vectors  and  White  Noise  Case 

In  addition  to  the  mean  vectors  being  equal  to  zero  on  the  two  hypotheses,  let 
the  noise  be  white  with  power  Nq.  Then 

E[n  n“]  =  No  I 

where  I  represents  the  identity  matrix.  Let  K,  denote  the  covariance  matrix  of  the 
signal, 

E[s  s«]  =  K.  . 

It  follows  that 
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Ko  =  E[f  f»|Ho]  =  E[n  n»]  =  I 


(5.9) 


and 

Kj  =  E[f  f^lHj]  =  E[(s  +  n)(s  +  n)^] 

=  E[s  s”]  +  E[n  n**]  =  K,  +  I. 

The  latter  result  uses  the  fact  that  the  signal  and  noise  are  uncorrelated  and 
both  have  zero  mean.  The  elements  of  the  vectors  s  and  n  ares^,  n^,  k=0,  1...,  N-1, 

respectively. 

Consider  again  the  inverse  relationship  between  and  K^, 

Qi  K,  =  1. 

It  was  shown  above  that 


K,  =  K.  .  No  1. 


Now  express  Qi  as 


Q  =  J_[i  -  H,]  (5.10) 

*  N  ^ 

where  is  an  N  x  N  matrix.  Substituting  for  Kj  and  in  the  first  equation  yields 


-  H,]  •  [K.  +  N„  n  =  I 
No 

which  simplifies  to 


K.  =  No  H,  +  K,.  (5.11) 

We  notice  that  this  equation  is  the  discrete-time  analog  to  the  integral  equation 
obtained  in  Sec.  2.4  for  the  continuous  case.  Since 
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we  have 


or  we  can  write 


Qo  =  =  [No  U-'  =  ^  I. 


<2i  =  Qo  ^  H, 


<Jo  Qi  =  ^  Hi 

No 

The  likelihood  ratio  test  can,  therefore,  be  expressed  also  in  terms  of  the  matrix  Hi: 


^R“H,  R  ^  Y 
N„  ‘  < 


(5.12) 


5.1.3  Comparison  with  an  Estimation  Problem 

Consider  the  minimum-mean-square-error  (MMSE)  estimation  problem  shown 

in  Fig.  5.1;  s,^g  =  f  denotes  the  minimum  mean  square  estimate  of  s  and 


f  =  s  +  n.  By  definition,  Sj 


s  and  f  must  be  orthogonal: 


E[(Smmse  ^'^”1 


(5.13) 


=  E[(HjS  +  Hjfi  -  s)(s  +  n)”]  =  0 


Assiuning  the  signal  and  noise  vectors  to  be  uncorrelated  and  zero  mean,  we  get 


Hi  K,  +  Hi  N„  I  -  K,  =  0 
K,  =  No  Hi  +  Hi  K,. 


(5.14) 


This  last  equation  is  the  same  as  the  one  foimd  previously.  We  conclude  that  Hi 
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characterizes  the  filter  which  generates  the  minimum  mean  square  estimate  of 


s  from  f. 


5.2  Ambiguity  Function  for  the  Discrete  Signal 


so  that 

s  =  6 

and 

K.  =  E[5  5«1  =  E  {l6p}  E^tt” 
=  2  oj  E,  T  t" 

=  E,  Y 


where  E^  denotes  average  energy 


We  now  introduce  the  matrix  Hi  from  Sec.  5.1.2  and  assume  Hi  to  be  of  the  form 


Hi  =  cY  y“ 

where  c  is  a  constant.  We  are  interested  in  finding  c  so  that  Hi  is  the  solution  of  the 
matrix  equation  derived  in  Sec.  5.1.2, 


K.  •  N,  H,  +  H,  K, 

which  is  the  disaete  form  of  the  integral  equation  in  Sec.  2.4.  Substituting  for  K,  and 
Hi  from  above  gives 


Ej  Y  Y*  =  Nq  c  Y  Y*  +  c  E  Y  Y^  Y  Y^.  (5.17) 

Notice  that  Y^  Y  is  a  scalar  and  represents  the  energy  of  the  vector  Y.  Assume  this 
energy  to  be  normalized  to  unity: 


We  have  then 


=  1. 


(5.18) 


E  Y  Y^  =  Nq  c  Y  Y^  +  c  E  Y  Y^ 
or 


E  =  N.  c+c  E , 

r  0  r' 


giving 


No-E, 


Therefore,  the  assumed  form  for  Hi  is  correct  and  Hi  is  given  by 


E  -  -H 

H  =  _ W  T  . 


The  likelihood  ratio  can  then  be  rewritten  as  follows: 

«,  =  —  R”  Hi  R 
^  N 

=  —  r“  T  R 

No  No  +  E, 

_  1  Ef  1 Rp 


1 

R”  Hi  R 

No 

1 

El 

No 

No  -  E, 

1 

El 

No 

No  ^E, 

1 

El 

No 

No  -  E, 

(5.19) 


(5.20) 


The  received  vector  R  can  be  written  as  the  sxun  of  the  incoming  actual  signal 


vector  s,  and  the  noise  vector  n: 


(5.21) 


R  =  s,  +  n  ' 

The  samples  which  constitute  assume  a  reference  delay  Xh-  Thus,  the  elements  of  s^ 
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=  6  f(kT^-cy““  k=0,l,...,  N-1  (5-22) 

where  t'  is  the  delay  discrepancy  -  Tg. 

For  the  parameter  estimation  problem,  the  likelihood  function  is  the  same  as 
the  likelihood  ratio.  At  first  we  assume  S  given,  so  the  likelihood  function  is 

.  p  +00 

Nq  Nq+Ej,  k=-<» 

C  +  00 

=  ^  At- I  5:  (b  Vo^y^'^r 

Nq  Nq+Ej  k=-<» 

=  -^1  E  (6  f(kT,-rO  e^"^*  +  V  (kT.)^ 

Nq  Nq+E,  k=-<» 

Here,  is  the  Doppler  shift  discrepancy  • 

Averaging  over  the  distribution  of  6  yields  the  likelihood  fxmction  as  a 
fimction  of  only  the  parameters  x'  and  cop; 

«R  =  /3(R|6)f6(6)d6  (5.24) 

where  fg(6)  is  the  density  function  for  8.  We  can  write 

y?  +  00 

«R  =  ^  ^  fl  ^  b  f(kT,-xV“^^-  ffc(b)  db 

No  No+E^  k=-<» 

+  (noise  x  noise  term)  +  (noise  x  signal  terms) 
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I  S  f(kT.-xOf’(kTy“°'^*P 

NoCNo^E^  k=-oo 


+  other  terms  involving  n^ 


(5.25) 


The  point  in  the  x',  plane  for  which  «r  is  maximum  represents  the  maximum 

likelihood  estimate  of  the  pair  (t^,  for  a  given  pair  of  values  (Xj,, 

The  continuous  time,  continuous  frequency  ambiguity  function  of  the  discrete 
signal  is  proportional  to  the  noiseless  part  of  the  likelihood  function  {r  and  is  given 


e.(V.o>;)  =  I  E  f(kr,-x^r(kT> 

k=-« 


(5.26) 


5.3  Relationships  Between  the  Continuous  Ambiguity  Function  of  the  Discrete  Signal 
and  the  Continuous  Ambiguity  Function  of  the  Continuous  Signal 

We  compare  now  the  AF  (ambiguity  function  )  of  the  discrete  signal,  obtained 
in  the  preceeding  Section,  Eq.  (5.26),  with  the  AF  of  the  continuous  signal  as 
presented  in  Sec.  2.3.,  Eq.  (2.48) 


0(x' =  |/7f(t-Tyf  (t)e^“*>‘  dti^ 


(5.27) 


Let  ^Tx',  0)^,)  denote  the  time-frequency  correlation  of  the  discrete  signal. 


+00  f 

k=-» 


(5.28) 
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so  that 


(5.29) 


Similarly,  let  be  the  time-frequency  correlation  of  the  continuous  signal  f. 


=  /7  fCt-T^lf  *(t)e^““‘  dt 


so  that 


=  I‘I>(t  ,Od) 


We  can  express  as  follows: 


+  00 

=  E  /'7(t-T^f’(t)e^^"^  6(t-kr>it 

k=  -oo  *" 

=  |"7(t-T')f*(t)e^“®*  E  «(t-kT^  dt. 

k=-<» 

Using  Poisson's  sum  formula,  the  sum  of  impulses  can  be  written  as  [7] 


+00  +00  -j2n-E-t 

E  fi(t-kTJ  =  -L  E  e 

k=-<»  Tg  p=-oo 


(5.30) 


This  results  in 


1  +* 

<l>d('f'.“D)  =  If  ^  *(t)e  dt 

r  _  __  i/  -• 


^  p=- 


+  00 

n=-oo  A* 


Finally, 


(5.31) 


114 


T,  p=-oo  • 


(5.33) 


Thus,  we  have  the  relationship  between  the  AF  of  the  discrete  signal  and  the 
continuous  signal. 
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6  DISCRETE-TIME  DISCRETE-FREQUENCY  AMBIGUITY 
FUNCTION  OF  A  TIME  SAMPLED  SIGNAL 


In  this  chapter  we  seek  an  expression  for  the  discrete-time  discrete-frequency 
ambiguity  function  of  a  sampled  signal.  We  begin  by  considering  a  signal  that  is 
sampled  both  in  time  and  frequency. 

6.1  Expression  for  a  Signal  Sampled  in  Both  Time  and  Frequency 

Consider  a  signal  s(t)  that  is  nonzero  only  on  the  interval  [0,Tf].  A  signal  §^(1) 

that  is  time  limited  to  interval  [0,7^]  can  be  obtained  from  s(t)  in  the  following 
mannen 


r  s(t)  OitiTj 

i  0  elsewhere 


where  <  Tf.  We  can  write  s,(t)  as 


where 


SeW  =  s(t).w^(t) 


(6.1) 


(6.2) 


Wd(t)  = 


1 


(6.3) 


[  0  elsewhere. 

A  periodic  version  of  s^ft)  is  now  created  by  replicating  s^(t)  every  1/X,  seconds 

where  the  spacing  between  the  replications  is  required  to  be  greater  than  or  equal  to 
Tj.  This  can  be  done  by  convolving  s,(t)  with  a  periodic  impulse  train,  as  shown 

below. 
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(6.4) 


Sp(t)  =  §,(!)♦  E  fiCt-y). 

m=-“> 


Notice  that 


§.(t)  = 


§(t) 


elsewhere. 


s  (t)  is  now  sampled  every  T  seconds  to  obtain  s^(t) .  Formally,  this  can  be 


expressed  as 


I  (t)  =  §  ft)  E  6(t-kT) 


=  E  s,(kT)  fi(t-kT). 
k=-«> 


Assxune  =  (N-l)T.  Note  that  s^CkT)  is  zero  for  k  <  0  and  k  >  N-1.  A  periodic 
version  of  I  (t)  is  obtained  by  replicating  i,(t)  every  ^  seconds.  This  results  in 


§p(t)  =  MO*  ^ 

m=-i» 


m=-<» 


+00  N-1 


=  E  E  s,(kT)  5(t-kT-f). 
m=-<»  k=0 


Note  that  I  (t)  is  discrete  in  time.  Because  it  is  also  periodic  in  time,  its  Fotuier 
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transform  is  discrete  in  frequency.  It  is  in  this  sense  that  ip(t)  represents  a  signal 

which  is  sampled  in  both  time  and  frequency.  Observe  that  the  m=0  term  in  the 
expression  for  ip(t)  is  equal  to  l^(t).  Since  s,(t)  is  equal  to  Sp(t)  for  t  €[0,  (N-l)T]Jp(t) 
can  also  be  written  as 

+00  N-1 

ip(t)  =  E  E  §00)  a(t-kT-^^). 

xn=-oo  k=0  ^ 

The  signals  arising  in  the  transformation  from  s(t)  to  lp(t)  are  illustrated  in  Fig.  6.1. 

Assume  that  the  spectrum  of  s^ft)  is  approximately  bandlimited  to  the  interval 

as  shown  in  Fig.  6.2.  Therefore,  to  prevent  aliasing  in  frequency,  the  sampling 
interval  T  should  be  chosen  such  that 

-  i  2B  (6.9) 

T 

(i.e.  the  sampling  rate  should  be  greater  than  or  equal  to  the  Nyquist  rate).  In  order 
to  avoid  aliasing  in  time,  the  sampling  rate  in  frequency  should  be  chosen  such  that 

1  i  T.  +  6  (6.10) 

X  ^ 

where  e  is  a  small  non-negative  constant  added  to  prevent  aliasing  of  the  end  points. 
Since  T^  =  (N-l)T,  note  that 

j-  i  (N-l)T+e  (6.11) 

or  equivalently 


9 


Se(t) 


-BOB 


Figure  6.2  -  Sketch  of  s^(t)  and  its  spectnun. 
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-  ^  NT. 
X 


(6.12) 


6.2  Discrete-time  Discrete-frequency  Ambiguity  Function 

The  continuous  time  continuous  frequency  correlation  hmction,  whose 
magnitude  squared  is  equal  to  the  AF,  is  given  by 


(6.13) 


where  the  subscript  f  denotes  the  fact  that  the  time-frequency  correlation  is  with 


respect  to  the  signal  f(  ).  Let 


g(t)  =  f(t-Tje^’'^‘. 


(6.14) 


g(t)  is  the  desired  target  return  at  the  input  of  the  radar  receiver  with  actual  time 

delay  x,  and  actual  doppler  frequency  f,.  f(t-T)e^’'*'  is  the  reference  signal  used  by 

the  matched  filter  in  the  receiver  with  reference  time  delay  x  and  reference  Doppler 
frequency  f.  For  convenience,  the  subscript  H  is  not  used  in  this  analysis.  As  a 

function  of  g(t),  we  can  write 


<|)|(x4)  =  /  g(t)f*(t-T)e-^"^ 


(6.15) 


This  choice  for  representing  <|)j(',0  is  meaningful  in  the  sense  that  it  involves  the 
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correlation  of  the  desired  target  return  with  the  reference  signal  to  obtain  the  time 
frequency  correlation  <|)j(-,  ). 

6.3  Sampling  the  signal  in  time  and  the  signal  spectnun  in  frequency 

Now  assiune  that  discrete  samples  are  present  at  the  input  of  the  receiver 
instead  of  a  continuous  waveform.  Let  the  discrete  samples  be  obtained  by  sampling 
a  continuous  waveform  both  in  time  and  frequency.  Following  the  same  procedure 

as  developed  in  the  previous  section,  the  versions  of  |(t)  and  f(t)  that  are  sampled  in 

both  time  and  frequency  are  expressed  as 

+00  +0O 

|p(t)  =  E  E  gJXT)  fi(t-kT-f) 

ni=-oo  k=-<» 

^  +00  +00 

yt)  =  £  E  f.{fr)  s{t-«T-f)  (6.16) 

Il=-00  {=—00 

where  the  sequences  g^CkT)  and  f,G^  have  been  limited  to  the  intervals  ['^T,  (N-l)T 
+  yNT],  Y  an  integer  and  [0,  (N-l)T],  respectively.  The  time  frequency  correlation 

A 

function  involving  |p(t)  and  fp(t)  is  given  by 
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+00 

— ee 

+«*  +<»  +00 

=  /  E  E  i.(kD  fi(t-kT-^) 

00  m=-<»  k=-<» 

+  00  +00 

•  E  E  f;(fiT)  5(t-T-«T-|) 

Il=-00  j=— 00 

Let  p=k-{  and  q=m-n.  Then,  «=k-p  and  n=m-q.  Note  that 

fi(t-kT-^)  8(t-T-n'-|) 

=  8(t-kT--)  Sa-T-kT+pT-^+f) 

=  6(t-kT-^)  8([t-kT-^]-[x-pT-f]) 

A  A  A 

=  fi(t-kr-^)  8(T-pT-|). 

Therefore, 


(6.17) 


(6.18) 


+00  +00  +  00  +  00 

=  r  E  •  r  E  8.(kT)f;(kT-pT) 

'  ni=-<»  k=-<»  q=-<»  p=-«> 

+  00 

f  5(t-pT-^)  8(t-kT-f) 

—  00 


+  00  +  00  +00  +00 

=  E  E  E  E  g,(kDf;(kr-pT). 

ni=-oo  k=-<»  q=-<»  p=-oo 


8(i:-pT-f) 


-PTlflkT^B] 

e 


Poisson's  sum  formula  states  that 


(6.19) 
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+  00 


+00 


(6.20) 


X)  e  ’■=  X  E  5(f-mA). 


m=-'» 


ni=-«> 


Thus, 


+00  +00  +  00  +00 

(x,f)  =  E  E  E  E  t,(kT)  fJ'CkT-pT) 

m=-<»  k=-<»  q=-oo  p=--oo 

«(t:-pT-|)  X  8(f-mA.)e-i2«»r. 


(6.21) 


Assume  X  aud  T  are  chosen  such  that  =  L  where  L  is  an  integer  greater  than  (N-1). 


Also,  let  V  =  =  p+qL.  Note  that  u  is  also  an  integer.  As  q  ranges  from  -<»  to  +«», 


u  also  varies  between  -o<>  and  +c«.  Introducing  the  change  of  index,p  =  v 


can  be  expressed  as 


+00  +00  +00  +00 

<l)j^(T,f)  =  X  E  E  EE  g,(kT)  f;(iT-vT+^) 

ni=-oo  k=-oo  q=-oov=-oo 

8(t-vT)  fi(f-mX) 


(6.22) 


By  defining  the  discrete-time  discrete-frequency  correlation  function  of  fp(t)  as 


+00  +00 


ki\>^)=x  E  E  g,(kT)f;(kr-vT+-a) 

^  k=-«q=-oo  ^ 


(6.23) 


the  expression  for  ("^4)  can  be  written  as 
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hOO  +00 


<K(x,f)=  ^  5^  ^(v,m)  8(x-vT)  8(f-mA,). 

'  —  — .AA  Tn~— AA 


(6.24) 


V  =  -<»  IIl=-oo 


Note  that  <|)j  (x^)  is  periodic  both  in  time  and  frequency. 


Sampling  only  the  signal  in  time 


From  the  expression  for  gp(t),  observe  that  the  m=0  term  results  in 


lp(t)|  g,(kT)5(t-iaO  =  i(t). 


(6.25) 


A 

Also,  from  the  expression  for  f  (t),  observe  that  the  n=0  term  results  in 


fft)  =  E  f,(«D6(t-«T)  =  f,(t). 
I  “=0  {=-00 


(6.26) 


Since  q=m-n,  m=n=0  implies  q=0.  Hence,  when  q=m=n=0. 


(6.27) 


By  straight  forward  substitution  of  the  expressions  for  g^(t)  and  fg(t-x),  we  obtain 


4>,(xd)  =  f  E  g,(kT)8(x-kD  E  f,(fD  5(t-x-«D  e-J2»«dt 

•  J  1 -  {=-00 


(6.28) 


_«  k=-oo 
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Let  v=k-L  Then  {=k-v.  Observe  that 

«(t-kT)  8(t-x-{T)  =  6(t-kl0  8(t-T-kT+vT) 

=  8(t-kT)  8[(t-kT)  -  (t-vT)] 

=  8(t-kT)  8(x-vT).  (6.29) 

It  follows  that 

00  00 

<|>|(x.f)  =  E  E  |^(kT)  f/CkT-vT) 

*  k=-<»  v  =  -<» 

00 

f  a(t-kT)  8(t-vT) 

—  00 

00  00 

=  t  E  I.TO  f;(KT-vT)  «(t-vT)  .  (5.30) 

k=-oo  v  =  -<“ 

Note  that  4^  (x^  is  sampled  in  the  delay  variable  x  since  the  correlation  function  is 

nonzero  only  for  x  =  vT  where  v  is  an  integer. 

Although  <|)f(x^  is  a  continuous  function  of  the  frequency  variable  f,  in 

practice  the  correlation  fxmction  will  be  evaluated  only  at  frequencies  spaced  X,  Hertz 
apart.  In  effect,  <|>f(x,t)  is  sampled  in  frequency  as  well.  Analytically  this  is 

accomplished  by  multiplying  <(>|  (x,f)  with 
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00 


(6.31) 


Consequentiy,  the  discrete-time  discrete-frequency  correlation  function  for  the  Ume 
sampled  signal  f^(t)  becomes 
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^  ^  f;(kT-vT)  e 

k=-<» 


-j2,t(^-gkr 


(6.36) 


This  is  the  form  we  will  use  for  the  discrete-time  discrete-frequency  correlation 
ftmction  of  the  truncated,  time  sampled  signal,  f^(t).  By  definition,  the  corresponding 
discrete-time  discrete-frequency  ambiguity  function  is 


6j^(v,m)  =  |^^(v,m)p  . 

Thus,  we  have  the  desired  expression  for  the  time  sampled  signal. 


(6.37) 
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7  LINK  BETWEEN  THE  DISCRETE-TIME  DISCRETE-FREQUENCY 
AMBIGUITY  FUNCTION  AND  THE  CONTINUOUS-TIME 
CONTINUOUS-FREQUENCY  AMBIGUITY  FUNCTION 

In  this  chapter,  we  are  interested  in  finding  the  relationship  between  the 
discrete-time  discrete-frequency  correlation  function  ^  (v,  m)  and  the  continuous-time 

continuous-frequency  correlation  function  (t,  f). 

7.1  Sampling  the  Signal  in  Time  and  the  Signal  Spectrum  in  Frequency 

By  definition,  the  signal  f^ft)  is  nonzero  only  over  the  interval  [0,  (N-l)T],  The 

periodic  signal  f^ft)  is  formed  by  replicating  f^(t)  every  seconds.  Specifically, 


*  co°*2(t) 

X 

where 

+  00 

comb^ft)  =  E  6(t-Y)- 

X 

fji(t)  is  obtained  by  sampling  f^ft)  every  T  seconds.  In  particular. 


fp(t)  =  yt)  •  combat) 

where 

+  0O 

combat)  =  E  5(t-kT). 
k=-<» 


(7.1) 


(7.2) 


(7.3) 


(7.4) 


129 


By  definition,  the  time-frequency  correlation  function  between  the  signals 


f(t)  and  |(t)  is 


When  |(t)  =  f(t-xj  4)jj(t4)  is  only  a  hmction  of  f(t).  Hence,  we  replace  the 


subscript  by  f  and  write 


<l>?(^4)  =  / g(t)  r(t-T)  e-J2”ftdL 


Consider  now  the  time-frequency  correlation  function  given  by 


=  /  ge(t)  fe(t-T) 


The  following  properties  hold.  These  are  derived  in  Appendix  A. 


1.  =  f <|>j^(p,x)<|) 

coDob  ifComb  1  (T-p,x)  dp 


Gooib  iiOonsb  I  (T-P,X)  dp 
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where 


gp(t)  =  ge(t)*  combi(t) 


(7.10) 


lp(t)  =  gp(t)  .  comb^Ct). 


(7.11) 


We  see,  therefore,  that 


^(tS)  =  //  tj/M)  dP"** 


(7.12) 


where 


V=-oo  f=-00 


+00  +00 


=  E  L  f  6(t-kT)  6(t-T-«T) 

k=-oo  {=-00 


(7.13) 
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+  00 


+00 


+  00 


^  «»-x)  ^  S(t-Ttp-i)e^*”dt 

_«,  {=-*  ]£=-«> 


X  X 


+  00  +00 


+00 


=  E  E  /  5(t-|)  5[t-(T-p)-^] 

{=-00  k=-<»  ^  * 


(7.14) 


Let  m=k-j.  Then  ®=k-m  and 


+00  +00 


^  ^  /  «(t-kT)8(t-T-kT.mT)e-»’»-«yt 

k=-00  m=-<»  _ao 


+  00  +00 


+00 


=  E  E  I  5(t-kT)  a(T-mT) 
k=-<»  m=-«>  _oo 


+  00  +00 

=  E  E  a(x-mT) 

k=-<»  m=-«> 


(7.15) 


Using  the  Poisson  siun  formula^. 


1  tr  -jZnBT 

E  a(T-mT)  =  -  E  e  ^  . 


m=-<» 


T  01= -00 


(7.16) 


Therefore,  the  expression  for  <{)  becomes 


+  00  +00 


(x,f-x)=l  E  E 


ccod>pconib].'  ■  '  rp 

^  k=- 


00  10=-OO 


(7.17) 


Following  a  similar  procedure  on  the  expression  for  ^  eonibi(^“P»*)  obtain 


X  A 


+00  +00 


-j2Rxi 


{=-00  n=-oo 


(7.18) 


X  X 
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Substituting  the  simplified  expressions  for  ♦^nbi.eombj^C'c-p^)  and<|>,^^(T^-x) 

k  k 

into  that  for  we  get 


Hence, 


(7.19) 


(7.20) 


(7.21) 
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=  -rr 


E  E  e"^8(x-p-“). 

f  =  — 00  11=  -00  ^ 


0®  ‘IB 

E  E  8(f-x--)  dpdx 

k=-oo  m=-«> 

00  00  OO  00  a 

_  Y,  E  e 

{=-00  n=-<»  k^-oo  in=-'» 


If  4>t(M  6(t-p-y)  8(f-x-|)  dpdx 


00  00  00  00  a 

—  E  E  E  E  e 


{=-00  n=-<»  k=- 


00  ni=-'» 


n  .  k,  -J2n(f-i)l 


(7.22) 


Recall  that  X  is  chosen  such  that  X,  =  1/LT,  where  L  is  an  integer  greater  than  N-1. 
Then 


J2n-«LT 

^  e  ^  =  1. 


(7.23) 


The  expression  for  <t>j  (x4)  is  simplified  further  by  again  recognizing  from  Poisson's 


sum  formula  that 


00  a  00 

E  =  T  E  5(T-mT) 


(7.24) 


This  implies  that  L  =  N.  Observe  that  this  value  for  L  is  greater  than  (N-1),  as 
required  to  prevent  aliasing  in  the  time  domain.  Hence, 
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+00  +00 


<|)f(T^  =  I  ^  E  hM  fiCx-mlO  ^(f~) 

p  T  ^  ^  p  NT 


(7.29) 


m=-oo  f=-ao 


where 


+  00  +00 

&(m,t)  =  ^  s  <|),  ([in-nN]T,  («-kN]-i-). 

.Ml  j^=_oo  lj^=-oo  * 


(7.30) 


^(m,0,  the  discrete-time  discrete-frequency  correlation  fimction  of  the  discrete 


periodic  signal  fp,  is  a  repetition  every  j  =  NT  seconds  in  time  and  ^  Hertz  in 

frequency  of  the  continuous-time  continuous-frequency  correlation  function  for  the 
continuous  time  limited  signal  f^(t). 

7.2  Sampling  the  Signal  in  Time  and  the  Correlation  Function  in  Frequency 

Equation  (7.30)  expresses  the  relationship  between  the  continuous-time 
continuous-frequency  autocorrelation  function  4>f  (',•)  and  the  discrete-time  discrete- 

frequency  autocorrelation  function  ^  (*,•).  Recall  that  the  signal  ^(t)  is  sampled  in 
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time  every  T  seconds  while  the  spectrum  of  yt)  is  sampled  in  frequency  every  X 
Hertz  where  A,  =  — . 

NT 

In  a  similar  manner,  we  now  find  the  expression  of  the  discrete-time  discrete- 
frequency  autocorrelation  function  ^  (%*)  of  the  signal  fjx)  which  is  sampled  only  in 


time  every  T  seconds.  Consider  the  time-limited  signal  yt)  is  obtained  by 


sampling  f^(t)  every  T  seconds.  In  particular. 


yt)  =  .  comb^ft) 


(7.31) 


where 


+00 


comb.r(t)  =  E  6(t-kT). 

k=-'» 


(7.32) 


Using  properties  developed  in  Appendix  A, 


+  00 


-00 


(7.33) 


where,  as  shown  in  the  previous  section. 
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+00  +00  „ 


(7.34) 


Tk=- 


00  in=-oo 


Substituting  the  expression  for  _ k^(t^-x)  into  that  for  <|)j(t^,  we  get 


+  00  +00 


=  f  ^  E  e 

*  T  k=-oo  m=-oo 


(7.35) 


From  Poisson's  sum  formula,  note  that 


•rw  too  - 

E  =  I  E  fi(f-x-^). 

k=-oo  T  ]j-=_oe  T 


(7.36) 


Hence, 


+  00  +00 


f(T,0  =  f  <l>f(T,x)  ^  E  6(f-x-^)dx 

_oo  A  k=-oo  m=-<»  ^ 


+  00  +00 


1  V  j.  f  f  K  -^”r 

—  L  E  (|)^W--)e 

1  k=-oo  m=-oo  ^ 


(7.37) 


This  expression  can  be  simplified  fxuther  by  recognizing  from  Poisson's  sum  formula 


E  e  ^  =T  E  fi(T-mT). 


(7.38) 


m=-'» 


m=-'» 


It  follows  that 
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(7.39) 


<|)|(x^  =  4  E  E  5(x-inT) 

*  T  k=-«.  m=-«  T 

+00  +00  - 

=  -  E  E  <|>j^(mT^--)  fiC-r-mT). 

T  k=-<»  m=-<»  ^ 

Since  f^(t)  is  time  limited  to  the  interval  [0,  TJ,  is  time  limited  in  the 

argument  x  to  the  interval  [t.-Tj,  x.+Tj],  as  shown  in  Appendix  B.  Observe  in  the 
latter  expression  for  <|)j  (t^E)  that  the  argument  x  takes  on  only  the  discrete  values  mT 

where  m  is  an  integer.  Due  to  the  finite  time  extent  of  m  needs  to  range  only 

from  y*(N-1)  to  y+(N-1),  where  y  is  the  integer  value  closest  to  in  order  to  cover  that 

portion  of  the  x-axis  for  which  is  nonzero.  Thus, 


,  00  n=Y+(N-l)  , 

4H(T,f)  =  -  E  E  <|>f(mT,f-^)  5(T-mT). 

*  T  k=-oo  m=Y-(N-l)  ’  ^ 


(7.40) 


Because  f^(t)  is  sampled  in  time  every  T  seconds,  we  see  that  4>f^(T,f)  is  sampled  in  x 


every  T  seconds.  In  other  words,  is  a  disaete-time  continuous-frequency 

autocorrelation  function. 

In  practice,  ^  (x^f)  will  be  evaluated  only  at  discrete  points  in  the  x-f  plane. 
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Consequently,  we  now  sample  <|>|  (T,f)  in  frequency.  Analytically,  this  is  accomplished 
by  multiplying  ^  (t;!)  with 


=  E  <t>j^(mT,aC-^);  -<»i{^8o,  Y-(N-l)imiY  +  (N-l).  (7.43) 

k=oo  *  T 


then 


Y+(N-1)  » 

<|^(T,f)  comb^(f)  =  —  E  E  fi(ir-mT)  (7.44) 

*  Tni=Y-(N-l)  e=-«  * 


Recognizing  that 


comb^(f)  =  K  E  5(f-«K), 


(7.45) 


it  follows  from  the  above  expression  that 


Y^(N-1)  ^ 

<jH(-c,f)  =  -  E  6(x-mT) 

^  m=Y“(N-l) 


(7.46) 


Note  that  the  k=0  term  of  ^  (m,C)  is  given  by 


^/m,C)|^^  =K  <t)j^(mT,{K);  -<»  ^  Y-(N-l)imiY  +  (N-l)- 


(7.47) 


Therefore,  the  k=0  term  is  nothing  more  than  a  constant  times  the  continuous-time 
continuous-frequency  autocorrelation  function  of  the  time  limited  continuous  signal 

f^(t)  evaluated  at  the  points  (mT,  flC)  in  Ore  t-f  plane  where  m  and  I  are  integers  such 
that  Y-(N-l)i:ms:Y+(N-l)  and  Also,  note  that  is  the  replication  in 

frequency  of  <|>j  (mT,«K)  every  ^  Hertz.  As  was  the  case  in  the  previous  section,  the 

DFT  of  N  time  samples  spaced  every  T  seconds  results  in  a  transform  whose 
frequency  components  are  spaced  every  ^  Hertz.  Consequently,  K  is  chosen  such 


K  =  — . 
NT 


(7.48) 
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Then  ^(m,0  becomes 


1  * 

^  <|)j/mT.y-kN]^);Y-(N-l)^m^Y+(N-l). 

IN  1  —00  * 

From  the  previous  section,  the  discrete-time  discrete-frequency  autocorrelation 

A  ^ 

frmction  of  the  signal  f^ft),  which  was  generated  by  sampling  both  f^(t)  in  time  and 
its  spectnun  in  frequency,  was  defined  to  be 


m.«)  = 


(7.50) 


It  is  seen  that  ^  (m,0  can  be  obtained  from  ^  (m,Q  by  letting  n=0  in  the  expression 


for  ^(m,Q.  In  a  similar  manner,  it  can  be  shown  that  the  discrete-time  discrete- 


frequency  autocorrelation  function  for  fit),  when  only  the  spectnun  of  f  (t)  is 


sampled  in  frequency,  can  be  obtained  from  ^(m,{)  by  letting  k=0  in  the  expression 


for  ^^(m,0,. 


7.2  Aliasing  considerations 

It  is  shown  in  Appendix  B  that  (T,f)  has  both  twice  the  time  dxuration  of  the 


signal  f^(t)  and  twice  the  frequency  dmation  of  the  spectrum  of  f^(t).  The  expression 
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that  relates  to  <|>i  (mT,[{-kN]:^)  shows  that  4)|(mT,-— )  is  replicated  in 

h  NT  * 

frequency  every  —  Hertz.  This  is  illustrated  in  Figure  7.1  where  the  spectrum  of  f^ft) 

is  assumed  to  be  bandlimited  to  the  frequency  interval  [-B,  B].  In  order  to  avoid 
aliasing,  it  is  necessary  that 

1  ^  4B.  (7.51) 

T 

Therefore,  the  signal  f^(t)  must  be  sampled  in  time  at  a  rate  that  is  no  less  than  twice 

the  Nyquist  rate  if  aliasing  is  to  be  avoided  between  the  replications  of<l>-^(mT,^). 

When  aliasing  is  avoided,  we  have 

<^(m,l!)  =  —  <l>f(niT,:^) 

NT  •  NT 

where  y-CN-l)  <  m  <  y  +(N-1).  From  Appendix  B  it  is  apparent  that  the  primary 
interval  is  spanned  over  frequency  for  those  values  of  «  such  that 
p_(N-l)  i  t  i  p+(N-l)  where  the  value  of  p  is  the  closest  integer  smaller  or  equal 

to  —  =  f  NT  Since  the  discrete-time  discrete-frequency  ambiguity  function  is  the 

A  ‘ 

magnitude  squared  of  the  discrete-time  discrete-frequency  time-frequency  correlation 
function  it  follows  that 


143 


(7.53) 


is  periodic  in  frequency  with  the  primary  interval  extending  for 

Y-(N-1)  ^  m  ^  Y+(N"1)  p-(N-l)  i  fi  i  p+(N-l). 

The  effects  of  aliasing  for  the  complex  envelope 

f(t)  ,  ^  ^  smc(2wt)  (7.54) 

27CWt 

are  shown  in  Figs.  7.2(a)  through  (d).  Because  of  the  presence  of  the  sine  function,f(t) 
is  bandlimited  to  a  bandwidth  of  w  Hertz.  By  analysis,  the  continuous-time 
continuous-frequency  AF  of  f(t)  is 


f-f.  (2w-lf-fJ)  .  U(\\i\2  (7.55) 

e^T,f)  =  I  lect  ( — - - — —  sme  [(T-Tj(2w-lf-fJ)]|  . 

^  4w  2w 

This  is  plotted  in  Fig.  7.2(a)  for  w=l.  The  discrete-dme  discrete-frequency  AF  of  the 
sampled  signal  f(t),  using  a  sampling  frequency  equal  to  twice  the  Nyquist  rate,  is 


shown  in  Fig.  7.2(b).  Note  the  periodicity  of  the  AF  along  the  frequency  axis  and  the 
lack  of  aliasing.  The  AF  in  the  primary  interval  of  Fig.  7.2(b)  is  identical  to  the 
continuous-time  continuous-frequency  AF  in  Fig.  7.2(a).  When  the  sampling  rate  is 
dropped  to  1.8  and  1.5  times  the  Nyquist  rate,  noticeable  aliasing  occurs,  as  shown  in 
Figs.  7.2(c)  and  (d).  The  aliasing  becomes  more  severe  as  the  sampling  rate  is 
decreased.  A  sampling  rate  less  than  twice  the  Nyquist  rate  can  be  used  only  if  the 
aliasing  that  is  introduced  can  be  tolerated. 
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Fig.  7.2  (b)  -  Discrete  ambiguity  function:  2  Nyquist  rate  sampling 
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Fig.  7.2  (c)  -  Discrete  ambiguity  function:  1.8  Nyquist  rate  sampling 
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Fig.  7.2  (d)  -  Discrete  ambiguity  function;  1.5  Nyquist  rate  sampling 


8  RECONSTRUCTION  OF  THE  CONTINUOUS  AF  FROM  THE 
DISCRETE-TIME  DISCRETE-FREQUENCY  AF 


In  this  chapter  we  consider  the  reconstruction  of  the  continuous  ambiguity 
function  when  the  discrete-time  discrete-frequency  ambiguity  fimction  is  available. 
A  number  of  cases  are  considered. 

8.1  Signals  Sampled  in  Both  Time  and  Frequency. 

A 

Recall  that  f^(t)  is  the  periodic  signal  obtained  by  sampling  both  f,(t)  and  its 
spectrum  in  time  and  frequency,  respectively.  Define 


rect^ft)  =  * 

^  0  elsewhere. 


Then  the  sampled  version  of  f,(t)  (see  Figure  6.1)  is  given  by 


f,(t)  =  yt)  rect2_(t)  . 


yt)  can  be  recovered  from  yt)  by  low  pass  filtering  the  latter.  In  particular. 


f,(t)  =  f,(t)  ♦  h„(t) 

where 

hj,(t)  =  sincC^). 

It  follows  that 


(8.3) 


(8.4) 
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f,(t)  =  E  f,(nT)  smc( 


t-nT, 


Recognizing  that 


f, (t)  =  fp(t)  recti(t)  =  l^(t)  f^(t)  (8.6) 

g. (t)  =  ^(t)  recti(t-T^  =  ^(t)  gjCt)  (8.7) 

and  using  the  product  property  (A.2)  from  Appendix  A,  the  time-frequency  correlation 

ftmction  of  f,(t)  can  be  written  as 


=  /  ♦^j.(P/-x)dx 


Mtee  ip(t)  =  yt-Tj  ,  yt)  =  recti(t),  and  gj(t)  =  recti(t-y.  WifliBe^tog,(t), 


notice  that  recti(t-xj  is  nonzero  for  i  t  i  +  Y 


Also,  recognizing  that 
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f, (t)  =  f.(t)  *  h^(t)  =  lit)  *  f^it) 

g. (t)  =  g,(t)  ♦  h^it)  =  lit)  *  gj(t) 


(8.10) 


and  using  the  convolution  property  (A.1)  from  Appendix  A,  the  time-frequency 
correlation  fimction  of  f^(t)  can  be  expressed  as 


=  /  ♦y,(PJD  *^(T-pj>ap 


=  /  1^(P4>  ♦^^(T-pjt)dp  =  4^(r4>  (811) 


where  ^(t)  =  and  g3(t)  =  bjp(t)  e^**^  .  gj(t)  arises  because^ 


in  order  to  pass  xuidistorted  the  spectrum  of  g,(t),  a  rect  function  centered  at  f,  should 
be  used  to  filter  the  spectrum  of  g,(t).  This  is  equivalent  to  convolving  gj(t)  with 

h^(t)  e^***‘.  We  see,  therefore,  that 


(8.12) 


where 
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=  f  rectj_(t-Tj  rectj^(t-p) 


(8.13) 


Fig.  8.1  shows  the  plots  of  rectj^(t-T^  and  rectj_(t-p).  The  product  of  the  two 

X  X 


functions  is  non-zero  for  the  following  cases: 

Case  A  P  <  +  —  *“id  t  runs  from  p  to  +  Y 


Case  B  <  P  <  «nd  t  runs  from  to  P  +  y 


For  Case  A,  <  P  <  +  -i.  Then 


P 

g-j2«cr-.)t  I  T.*l 

"  -j2a(f-x)  I  9 

»  rr  «(r-*)K4> 

g-p»(f-l)p  _  g  * 

j2n(f-x) 


['^  -  i]  . ,  -  s] 

j2a(f-x) 


(8.14) 


For  Case  B,  x^-y  <  P  <  V  Then 
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-j2it(f-x) 


j2n(f-x) 


rp*^t  1 1  -j2«(f-*) 

,-»•«-«  [-T^ir]  e 


(8.15) 


j2ii(f-x) 


Noticing  that  P  >  for  Case  A  while  P  <  ‘c^  for  Case  B,  the  results  obtained  for 
^  (p^-x)  for  both  cases  A  and  B  can  be  combined  to  result  in 


♦^(P,f-»)  =  -  e 


-ja-M 


sin  I  2ic(f-x) 


■|M.I  1 


2  2X\ 


ir(f-x) 


(8.16) 


The  expression  of  “PjO  i®  given  by 


(8.17) 


It  is  readily  shown  [4]  that  ♦-,(t-P4)  can  also  be  expressed  as 


=  /  M’-iV  dv 


(8.18) 


where  Hi-(jv)  is  the  Fourier  transform  of  hjj(t).  Since  hu(t)  is  the  sine  function,  its 
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Fourier  transform  is  given  by 


=  T  recti  [T(f+^)]. 

T  QiT) 

The  expression  for  becomes,  therefore, 

4-OB 

=  T*  r  recti [T(v-f,+^)]  recti [T(v-f +-?-)]  ®  dv. 

^  T  PT)  j  (21^ 

Sketches  of  recti  [T(v-f^+^^)]  and  recti  [T(v-f+^^)]  are  shown  in  Fig. 

product  of  these  two  functions  is  nonzero  for  the  following  cases: 

Case  A  f  <  f  <  f-+—  and  v  runs  from  f-— i-  to  f  +-^ 

*  T  aT>  *  (?T) 

£?se  B  f  -1  <  f  <  f  and  v  runs  from  f,-—  to  f+— 

■  T  •  •  crn  ai) 

For  Case  A,  f  <  f  <  f  +1.  llien 

•  *  T 


(8.19) 

(8.20) 

8.2.  The 
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f  +-i_ 

•  (2T) 

dv 

f-JL 

(2T0 

j2x(x-p)  I  f-^ 
j2x(T-P) 

g  2  ®  1? 

j2x(T-P) 


(8.21) 


For  Case  B,  f  -  1  <  f  <  f  Then 

■  Y  A 


♦^(x-P^ 


f+— 

(21) 

J  T*  »  dv 


T2  I 

j2ll(T-P)  I 

j2*(T-» 


=  T2 


J2»(W.-|)  (^)  -J2«(f-r,*i)  (^) 

e _ ^  ^  -  e  ^  * 

j2x(x-P) 


(8.22) 


Notice  that  f  >  £»  for  Case  A  while  f  <  f,  for  Case  B.  The  results  obtained  for 
♦gi^(T-p^  for  both  Case  A  and  B  can  then  be  combined  to  result  in 


=  -  t»  . 


x(x-p) 


(8.23) 
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Substituting  the  expression  derived  for  and  -P4)  *at  ot^ix^ 


we  obtain 


=  //  T*  ♦f/M)  • 


sin  |2a(f-x) 


IP-^ml  1 


2  2AJ 


sm|2ic[lf-f.|-Y]^] 

*(x-P) 


^-J2,(f-.)[^.i]  (8.24) 

We  now  replace  the  expression  for  ^(P^),  as  derived  in  Section  7.1,  Eq.  (7.29). 


Also,  letting  JL  =  ^  becomes 


+a»  +«■  I 

♦,(t4)  =  T*//  ^  ^  »(P-mT) 

^  *  m  -  -OP  |=— o» 


sin|2x(f-x) 


1 


2  2X1 


lt(T-P) 


(8.25) 


Carrying  out  the  integrations,  we  have 
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i,m  =  f.(«)  •  M*) 

>  f.(t)  .  due 

+09 

=  E  f,(nT)  sine  .  (8.27) 

n=-«»  \  T  / 

Using  property  (A.1)  of  Appendix  A,  we  can  write  that 
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(8.28) 


where  ♦_»(t“P4)  was  derived  in  Section  8.1,  Eq.  (8.23).  It  is  given  by 


♦^(T-P.0  =  -T^  . 


2  sin  |f  f^l  ^  ^ 


n(T-p) 


(8.29) 


Therefore, 


♦^(t4)  =  /  -  e 


IC(T-P) 


dp  . 


(8.30) 


Recall  from  Section  7.2,  Eq.  (7.46)  that 


1  Y+(N-1)  , 

♦f(P4)  =4  ^  »0-™T) 

*  ^  m=Y-(N-l) 

where  y  is  the  closest  integer  to  x/T.  Substituting  the  expression  for 
♦j  (p^  into  that  of  the  latter  becomes 


(8.31) 


frequency  correlation  function  of  the  time  sampled  signal 

In  a  previous  section  it  was  shown  that  aliasing  in  the  discrete-time  discrete- 
frequency  AF  can  be  avoided  by  sampling  f^(t)  at  a  rate  greater  than  or  equal  to  twice 

the  Nyquist  rate.  When  the  Doppler  frequency,  f^,  is  known,  the  Nyquist  rate 
suffices,  as  explained  below.  The  Doppler  shifted  spectrum  of  f,(t)  is  sketched  in  Fig. 

8.3  for  f,  >  0.  Since  f,  is  known,  it  is  possible  to  shift  the  spectrum  to  the  origin,  as 
shown  in  Fig.  8.4,  by  appropriately  mixing  the  received  signal  with  a  sinusoidal 
oscillator  tuned  to  f,.  For  the  spectrum  of  Fig.  8.3  the  Nyquist  rate  is  given  by 

fs  =  2B.  (8.33) 

Samples  obtained  employing  the  Nyquist  rate  can  be  used  to  interpolate  to  any  other 

values  for  f^(t).  To  avoid  aliasing  of  the  discrete-time  discrete-frequency  AF,  recall 
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Spectrum  of  fg  ( t) 


Fig.  8.3  Doppler  shifted  spectrum  of  f,(t)  assuming  f,  >  0. 
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Shifted  Spectrum  of  fg  ( t) 


Fig.  8.4  Spectrum  of  Fig.  8.3  shifted  to  the  origin 
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that  the  sampling  rate  should  be  at  least  as  large  as 

2f.  =  4B.  (8.34) 

However,  using  interpolation,  the  4B  samples  per  unit  time  can  be  generated  from  the 
2B  samples  per  unit  time.  Hence,  if  the  Doppler  frequency  is  known,  aliasing  can  be 
avoided  even  when  the  sampling  rate  is  as  small  as  the  Nyquist  rate. 

When  the  Doppler  frequency,  is  imknown  but  its  maximum  value,  ^  is 

known,  alaising  can  be  avoided  by  sampling  at 

f,  =  2(B+f,_)  (8.35) 

where(B+f,  ^)  is  the  largest  possible  frequency  for  the  Doppler  shifted  spectrum  of 
f^(t)  This  is  true,  because,  as  was  done  for  the  known  Doppler  case,  the  2(B+f,  ^) 

samples  per  unit  time  can  be  used  to  generate  the  4B  samples  per  unit  time  required 
to  avoid  aliasing. 

When  the  Doppler  shift  is  completely  unknown,  it  is  guaranteed  that  no 
aliasing  will  occur  only  if  the  sampling  rate  is  greater  than  or  equal  to  4B,  which  is 
twice  the  Nyquist  rate. 
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9  SIGNAL  REPRESENTATION  FOR 
THE  NON-COOPERATIVE  CASE. 


9.1  Introduction 

From  a  signal  processing  perspective,  non-cooperative  bistatic  radar  differs  from 
the  cooperative  case  with  dedicated  transmitter  in  two  major  respects: 

i)  The  transmitted  signal  is  not  designed  specifically  for  bistatic  radar 
applications,  and 

ii)  A  significant  discrepancy  can  exist  between  the  actual  transmitted  signal  and 
the  signal  which  at  the  receiver  is  assiuned  to  have  been  transmitted. 

The  implications  of  (i)  are  a  generally  reduced  performance  compared  to  the 
cooperative  case  with  dedicated  transmitter,  and  the  possible  need  for  special 
processing  which  takes  into  account  the  nature  of  the  transmitted  signal  -  e.g.,  TV 
broadcast  signals,  FM  broadcast  signals.  These  issues  require  more  extensive 
investigation  and  will  not  be  addressed  here. 

The  implications  of  (ii)  are  also  a  reduction  in  detection  performance  compared 
to  the  cooperative  case  with  dedicated  transmitter.  A  first  step  toward  a  thorough 
evaluation  of  this  effect  is  a  suitable  signal  representation  for  the  non-cooperative 
case.  This  is  considered  in  Section  9.2. 

9.2  Signal  Representation 

We  consider  an  actual  transmitted  signal  with  complex  envelope  as  in  (2.12), 

f(t)  =  A(t) 

resulting  in  a  received  signal  as  in  (2.22),  asstuning  a  single  point  target: 
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SrW  =  4  Re  {  ^6f(t-c.)  e^"“‘  e'"-'}. 


(9.2) 


An  incorrect  representation,  at  the  receiver,  of  the  transmitted  signal  can 
involve  errors  in 

a)  the  assumed  transmitted  carrier  frequency 

b)  the  epoch,  i.e.,  the  exact  timing  of  the  transmitted  signal 

c)  the  complex  envelope 

We  see  from  (9.2)  that  an  error  in  co^,  which  will  be  denoted  cOg,  adds  to  or 
subtracts  from  cod*-  Similarly,  an  error  in  epoch,  denoted  Xg,  adds  to  or  subtracts  from 
Xg.  Specifically,  suppose  that  the  reference  waveform  available  at  the  receiver  is 
delayed  by  Xg  compared  to  the  actual  transmitted  signal,  and  up-shifted  in  earner 
angular  frequency  by  cOg,  giving 

f^t)  =  A(t-T^)  =  f(t-T^ 

Expressing  (9.2)  in  terms  of  f^^f(t)  as  given  in  (9.3),  we  have 

Sg(t)  =  Re  {/e;6 

where  any  complex  constants  are  absorbed  into  b. 

Consider  now  an  error  in  the  complex  envelope.  Let  the  transmitted  complex 
envelope,  as  assumed  at  the  receiver,  be 
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f(t)  +  fe(t)  =  A(t)  +  A^(t)  (9.5) 

The  term  is  assumed  to  be  devoid  of  a  component  proportional  to  f ;  that  is, 

/  f(t)  =  0. 

Combining  (9.5)  with  (9.3)  gives  for  the  reference  waveform  at  the  receiver 

In  terms  of  this  reference  waveform,  the  expression  for  the  received  signal  (9.2) 
becomes 

Sr®  =  Re  { 

Next  we  consider  how  the  various  errors  reflected  in  (9.7)  affect  the  expression 
for  the  ambiguity  function.  We  use  the  sjonmetrical  form  of  the  ambiguity  function 
given  in  (2.51), 

=  I /■  f ‘(t+^)  e'"”"  dtp,  <9® 

where  and  (Od  “  “d  ~^d  •  In  order  to  express  the  ambiguity  fxmction 

as  a  function  of  Rr  and  the  receiver-to-target  range  and  the  target  velocity 
component  along  the  bistatic  bisector,  the  substitutions  derived  in  Chapter  4  have  to 
be  used.  This  is  not  done  here  in  order  to  simplify  the  discussion. 

Now  let  the  received  signal  (9.7)  get  processed  with  hypothesized  delay  and 
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Doppler  Tjj  and  COp  ,  respectively.  Thus,  let 

H 


and  (Od  =  Qd -w^-tOn 


// 


This  leads  to  the  following  expression  for  the  ambiguity  fimction: 


=  I  /■  (t-^)  dt|» 


(9.9) 


(9.10) 


The  infinite  limits  on  the  integral  have  been  carried  along  but  may  not  be  applicable, 
depending  on  the  nature  of  the  transmitted  waveform  (see  the  discussion  leading  up 
to  (2.42)).  We  see  that  the  time-frequency  correlation  function  in  this  case  consists  of 
the  time-frequency  correlation  of  the  reference  signal,  and  an  error  term: 


<|.(T",e.>^  =  /•  Ut-^)  dt 


(9.11) 


Equations  (9.10)  and  (9.11)  express  the  ambiguity  fimction  and  time-frequency 
correlation  function  from  the  point  of  view  of  the  receiver,  that  is,  based  on  An 


alternative  formulation,  in  terms  of  the  transmitted  signal  f ,  may  be  more  useful  for 
system  design  and  analysis  purposes.  This  results  in  the  expression: 
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(9.12) 


e(T",<o£)  =  I  +  f;(t+|)]  e'"^’ dtp 

4.(x"o>^  =  /■  f(t-:^)  r(t.|:)  dt 

J  -CP  2  2 

+  /“  fd-f)  e'*^'  dt. 


(9.13) 


We  observe  that  the  time-frequency  correlation  function  in  (9.11)  and  in  (9.13) 
consists  of  the  superposition  of  two  separate  functions.  Concentrating  on  the 


formulation  (9.13),  (Op)  is  the  sum  of  the  ordinary  time-frequency  correlation 


of  the  transmitted  signal,  <l)„(T^^,C0p),  and  an  error  term,  (|)g(T^^,(Op).  Both  these  terms 

incorporate  the  errors  in  epoch  and  in  frequency,  and  (0^.  Thus,  (9.13)  can  be 
written 

<|)(x^(0d)  =  <j)tt(x",(0D)  +  <j)^(T'',(0^. 

The  ambiguity  function  (9.12)  then  becomes 

e(x",o>^  =  |<|..(x",9>^  X  <|,,(x'>^P 

(9.15) 

=  0jj(x^(Od)  +  6^(x'',cod)  +  2Re  (j)^ 

where 
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e.(T",<0^  =  I  /  r(t+^)  e^“"‘dtp. 

—  00 

00 

=  I  /  f(f~)  f>i)  0“^'  *P- 
—  00 

and  we  have  the  bounds 

e,(T".0)")  +  0^(t",qS  -  s  0(t",(o^ 

i  0,(t",<o^  ^  0,(t",O)£)  +  2^, 

Qearly,  the  effect  of  an  error  in  complex  envelope  is  somewhat  unpredictable, 
and  this  type  of  error  needs  further  study.  For  instance,  the  nature  of  the  transmitted 
waveform  can  significantly  influence  the  magnitude  of  this  type  of  error. 

9.3  Discussion 

What  preliminary  conclusions  can  be  drawn  from  the  development  in  Sec.  9.2? 
A  general  form  of  the  ambiguity  function  for  non-cooperative  bistatic  radar  is 

given  in  (9.12),  referenced  to  the  transmitted  signal  f  and  incorporating  errors  in  the 
reference  waveform  used  by  the  receiver,  The  simplest  case  arises  when  the 

difference  between  and  f  is  solely  due  to  a  timing  error  and  an  error  in  carrier 

frequency,  In  that  case,  the  term  f*(t+Y)  within  the  bracket  is  zero,  and  6(t"  co^ 

has  the  same  shape  as  the  ordinary  ambiguity  ftmction  for  the  signal  which  is 
radiated  at  the  transmitter.  However,  the  function  suffers  erroneous  displacements 
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in  the  x  and  cDd  directions  by  amounts  and  respectively.  In  (9.12)  this  is 

evidenced  by  the  fact  that  the  variables  are  x"  and  (as  defined  in  (9.9)),  instead  of 

x'  and  ti)^.  The  manner  in  which  these  displacements  in  x  and  co^  affect  the  bistatic 

ambiguity  function  (expressed  as  a  function  of  Rr  and  V,)  will,  of  course,  depend  on 
the  bistatic  geometry  and  can  be  analyzed  using  the  principle  of  Chapter  4. 

If  the  discrepancy  between  and  f  cannot  be  expressed  in  terms  of  errors  in 

timing  and/or  frequency,  we  have  an  "error  in  the  complex  envelope".  This  case  is 
much  more  complicated.  The  time-frequency  correlation  function  (9.13)  now  has  an 

extra  term  due  to  the  error  component  f^,  which  enters  in  a  nonlinear  manner  into  the 
ambiguity  function. 

Generally  speaking,  this  case  will  give  rise  to  a  spreading  of  the  ambiguity 
function.  A  detailed  analysis  of  this  case  will  have  to  consider  specific  transmitted 
waveforms,  specific  bistatic  configxuations,  as  well  as  the  likelihood  of  various 
amounts  of  error  in  the  complex  envelope  for  a  given  type  of  transmitted  signal. 

Finally,  if  in  addition  to  an  error  in  the  complex  envelope  there  also  are  timing 
and  frequency  errors,  then  the  ambiguity  function  with  distortions  just  described  is 
subjected  also  to  the  displacements  in  the  x  and  cOp  directions,  as  in  the  first  case. 
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10  SUMMARY  AND  CONCLUSIONS 

This  investigation  was  primarily  devoted  to  the  ambiguity  fxmction  analysis  for 

bistatic  radar  systems.  The  major  contributions  of  this  effort  are  outlined  below. 

1.  In  most  of  the  existing  literature,  ambiguity  fxmction  is  simply  defined  without 
providing  any  relationship  with  fundamental  radar  signal  processing  problems. 
We  showed  the  manner  in  which  the  ambiguity  function  arises  naturally  when 
solving  the  radar  signal  detection  problem  and  the  associated  parameter 
estimation  problem.  The  derivation  of  ambiguity  function  based  on  the  basic 
principles  of  detection  and  estimation  theory  was  presented. 

2.  The  manner  in  which  the  bistatic  geometry  affects  the  ambiguity  function  was 
derived.  It  was  shown  that  delay  and  Doppler  are  no  longer  the  appropriate 
arguments  for  plotting  the  ambiguity  fxmction  due  to  the  geometry  and  the 
nonlinear  relationships.  Illustrative  examples  showed  the  unexpected 
behavior  of  the  ambiguity  function  plots. 

3.  For  sampled  radar  signals,  expressions  for  both  continuous-time  continuous- 
frequency  and  discrete-time  discrete-frequency  ambiguity  functions  were 
derived.  Relationships  between  them  were  established.  It  was  shown  that  a 
sampling  rate  of  greater  than  or  equal  to  twice  the  Nyquist  rate  is  needed  to 

avoid  aliasing. 

4.  A  signal  model  that  incorporates  various  xmcertainties  in  the  reference  signal 
at  the  receiver  for  the  non-cooperative  bistatic  radar  was  developed. 
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APPENDIX  A 


TWO  PROPERTIES  OF  THE  CORRELATION  FUNCTION 


In  this  appendix,  two  properties  of  the  time-frequency  correlation  function 
needed  for  the  analysis  of  the  main  text,  are  derived. 

Recall  that 

+  00 

-00 

Property  1  (convolution  of  two  time  signals) 

If  f(t)  =  fi(t)*f2(t)  and  g(t)  =  gi(t)*g2(t)»  where  the  asterisk  denotes  convolution, 

then 


(A.l) 


Proof 

The  left  hand  side  of  (A.1)  can  be  written  as 
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=  ///  iiW  fi‘(y)  g2(t-x)  ©■^“'‘dxdydt. 

-00 


The  right  hand  side  of  (A.1)  is. 


/  gi(u)  f,‘(u-A,)  e  ‘^"^  du 


{  gjCv)  ^*(v-T+;i)  e-J2’’^dv 


d;L 


fff  ii(u)  fi*(U“^)  ^2^^)  f2*('^“T+>.)  e dudvd^.. 


In  the  above  expression,  let  t=u+v  so  that  v=t'U.  The  above  triple  integral  then 
becomes 


00 

///  giC'i)  fc(t-u)  f2*(t-u-T+A,)  dudtdA. 

-00 

Now  let  x=u  and  y=u-X,  The  triple  integral  then  reduces  to 


00 

///  giW  fi*(y)  fcCt"*)  f2*(t-T-y)  e'^'^didydl 
—  00 


which  is  the  same  as  the  left  hand  side  of  the  expression  found  previously. 


product  of  two  time  signals) 


If  f(t)  =  ffi)  fjCt)  and  g(t)  =  |i(t)  g^it),  then 


(A2) 


Proof 


By  definition. 


-00 

This  represents  the  left  hand  side  of  (A.2).  Let's  now  expand  the  expression  on  the 
right  hand  side.  We  obtain 

00  «  * 

J  <|»j^  fj(x,x)  dx  =  I  f  giftj)  fi(trt)  e  ^  dtj 

_oo  -«>  -00 

-00 

'  Iff  ifio  fcw  *,  dt,  &, 


However, 


f  dx  =  6(t2-ti). 


Therefore,  the  right  hand  side  of  (A.2)  becomes 
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+  00 


5(tj-t,)  dt,  dtj 


+00 


=  /  8i(t,)  fed,)  e  '"*"*  dti 


which  is  the  same  as  the  expression  foimd  for  the  left  hand  side  of  (A.2). 
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APPENDIX  B 


TIME  AND  FREQUENCY  EXTENT  OF 
THE  CORRELATION  FUNCTION 

By  definition,  f^(t)  is  time  limited  to  the  interval  [0,  Tj]  where  Tj  =  (N-l)T  is 
the  time  duration  of  f,(t)  and  N  equals  the  number  of  time  samples  uniformly  spaced 

by  an  interval  T  over  the  signal  duration.  Assume  that  f,(t)  has  a  spectrum  that  is 

approximately  bandlimited  to  the  frequency  interval  [-B,  B]  where  B  is  defined  to  be 
the  signal  bandwidth  as  shown  in  Fig.  B-1.  To  avoid  aliasing  of  the  frequency 

spectrum,  the  sampling  rate  for  f^(t)  should  be  greater  than  or  equal  to  the  Nyquist 
rate.  Since  T  denotes  the  sampling  interval,  the  Nyquist  principle  requires  that 

-  i  2B. 

T 

We  now  consider  sampling  in  frequency  the  spectrum  of  f^ft).  Uniformly 
sampling  in  frequency  with  a  sampling  interval  equal  to  A,  Hertz  causes  f^(t)  to  be 
replicated  in  time  every  seconds.  To  avoid  aliasing  in  time,  it  is  necessary  that 

y  i  T^  =  (N-l)T+e 
where  e  is  a  small  positive  constant.  Equivalently, 
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fe(t) 


Spectrum  of  fg  ( t) 


Fig.  B-1  Sketch  of  f,(t)  and  its  spectrum. 
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-  i  NT. 
k 

Recall  that  the  continuous-time  continuous-frequency  autocorrelation  function 
of  f^(t)  is  given  by 

00 

=  /  g,(t)  f;(t-T)  dt 

—00 

where 

From  Schwartz's  inequality, 

00 

s  /  is.(t)iif;(t-T)i  dt 

-00 

00 

-00 

Sketches  of  |f,(t-T^|  and  lf,(t-T)|  are  shown  in  Fig.  B-2.  By  inspection  of  Fig.  B-2,  the 
integrand  in  the  bound  on  |<|>j^(t4)  |  is  nonzero  for  ^  t  only  when 

<:  T  <:  T,  +  Tj. 

Note  that  this  is  a  time  interval  of  duration  27^  centered  at  x,.  Because  of  the  time- 
limited  duration  of  it  is  seen  that  4>j(T,f)  is  also  time  limited  with  a  duration 

equal  to  twice  that  of  f^(t).  In  addition,  Schwartz's  inequality  reveals  that  has 
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its  maximum  value  when  x  =  x,  and  f  =  f,. 

Having  determined  the  duration  of  4^  (x^  along  the  x  axis,  we  now  examine 

its  extent  along  the  frequency  axis.  It  is  readily  shown  that  can  also  be 

expressed  as 

00 

=  f  G„(jv)  F;ij(v-f)]  6^2--  dv 
^00 

where  F^Qv]  is  the  Fourier  transform  of  f^(t)  and  G^ljv]  is  the  Fourier  transform  of  g^(t) 
given  by 

6.av]  -  F.B(v-y]  e"*""-. 

From  Schwartz's  inequality, 

+00 

|♦^(T4)|  i  /  |G.Uv]|  |F;[j(v-f)]|  dv 
-00 

+  00 

=  /  |F.B(v-f.)]|  |F.B(v-f)ll  dv  . 

-00 

Sketches  of  iF^|j(v -f,)]  [  and  lF,[j(v-f)l|  are  shown  in  Fig.  B-3. 

By  inspection  of  Fig.  B-3,  the  integrand  in  the  boimd  on  l<j>j^(T,f)|  is  nonzero 
for  aU  V  only  when 

f,  -  2B  <  f  <  f,  +  2B. 

Notice  that  this  is  a  frequency  interval  of  duration  4B  centered  at  f..  Because  the 
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0  f-  B  f  f+B 

Fig.  B-3  Sketches  of  |F,[j(v-g]  I  and  |F,|j(v-f)]|. 
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spectrum  of  f,(t)  is  confined  to  the  interval  it  is  seen  that  <|>|  (t^)  is  confined 

to  a  frequency  interval  twice  as  large.  As  mentioned  previously,  Schwartz's  inequality 
shows  clearly  that  <j)|  (t^  has  its  maximum  value  when  i  =  t,  and  f  =  f,. 

In  smnmary,  <|)}  (t^  is  centered  at  and  is  limited  both  in  time,  with  a 

duration  equal  to  twice  that  of  f^(t),  and  in  frequency,  with  a  frequency  extent  equal 

to  twice  that  of  the  spectrum  of  f^(t). 
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